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In  this  report  we  discuss  adaptive  phase  compensation  for  the] 
effects  of  turbulence  and  thermal  blooming.  Closed-loop 
systems  employing  both  outgoing-wave  and  return-wave  control 
algorithms  are  considered,  although  we  have  obtained  from 
quantitative  estimates  of  system  performance  only  for  the 
return  wave- type  of  Bystcm.  On  the  basis  of  results  obtained] 
from  a computer  simulation  of  such  systems,  W"  conclude  that 
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20.  AHSTRACT  (continued) 

[ adaptive  optical  Byatems  employing  a return-wave  control 
policy  will  yield  only  limited  compensation  for  the  effects 
of  cw  thermal  blooming  unless  the  target  has  glints  that  ^ry 
significantly  smaller  than  the  diffraction-limited  width  or 
the  beam.  We  explain  this  in  terms  of  an  isoplanatic  effect, 
in  which  the  phase  in  forma:,  ion  returned  from  the  glint  is 
garbled  by  virtue  of  its  finite  extent.  In  contract  tc  the 
cw  thermal  blooming  results,  wo  find  that  return-wave  sys- 
tems are  much  less  sensitive  to  glint  size  for  turbulence  and 
nultipulso  thermal  blooming  tompenaiiition  (70  to  'JO*  of  free- 
space  peak  11'radi.anoe  reatorid  typically  in  the  case  of 
turbulence  eompenr.it ion  and  n facto v of  two  tc  thve.?  increase 
in  target  i i rad i ance  in  the  r a&'e  c f.  i nui.tipulue  thermaJ.  bloom- 
ing). Analytical  and  compute*  si  tuj  1 1 tlon  ’ results  have  also 
been  obtained  for  the  depenchs;  ce  of  m irror  fit  on  the  number 
of  mirror  actuators.  > 
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SUMMARY 


In  this  report-,  we  discuss  adaptive  phase  compensation 
for  the  affects  of  turbulence  and  thermal  blooming.  Closed- 
loop  systems  employing  .both  outgoing-wave  and  return-wave 
control  algorithms  are  considered,  although  we  have  obtained 
quantitative  estimates  of  system  performance  only  for  the 
return  wave-type  of  system.  On  the  basis  of  results  obtained 
from  a computer  simulation  of  such  systems,  we  conclude  that 
adaptive  optical  systems  en, cloying  a return-wave  control 
policy  will  yield  only  limi -.ed  compensation  for  the  effects 
of  cw  thermal  blooming  unl-eis  the  target  has  gMnts  that  5, re 
significantly  smaller  than  fcho  diffraction-limited  width 
of  the  beain.  . V03  oKpl^in  chin  . 1 terms  of  an  in  5pJ.au  e\<: 

1 1 

effect,  ir.  \wllirl-'  the  phase  it  formation  returned  frc*  tl.'n  flint, 
is  garbled  lr>.y»  ’’irt.ie  cf  ir.11  t.V.n.lte  extent , In  contrail  1.  th*> 
cw  thermal  blooming  results,  we  find  that  return -wave  ay  it.  terns 
are  much  less  ue.naitivn  to  glint  ai:iv.  :'or  turbulence  and 
multipulse  thermal  blooming  compensation  (70  to  90%  of  free- 
space  peak  irradiance  restored  typically  in  the  case  of 
turbulence  compensation  and  a factor  of  two  to  three  increase 
in  target  irradiance  in  the  case  of  multipulse  thermal 
blooming) . Analytical  and  computer  simulation  results  have 
also  been  obtained  for  tne  dependence  of  mirror  fit  on  the 
number  of  mirror  actuators. 
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INTRODUCTION 
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The  principal  objective  of  moat  high  energy  laser 
systems  is  to  achieve  the  maximum  possible  irradiance  at  a, 
distant  target.  Under  ideal  conditions  the  maximum  irradiance 
if  determined  by  free-spac©  diffraction  effects.  In  practice, 
however , this  ideal  state  is  seldom  realized.  We  must  contend 
additionally  with  wave  distortions  introduced  by  the  lauer 
device,  by  the  optical  train  and  by  medium  inhomogeneitiua 
in  the  path  between  the  transmitter  aperture  and  the  target. 

In  each  case  the  * di  stortion  is  caused  larguly  by  path  length 
errors  that  introduce  phase  aberrations.  In  this  report, 
we  consider  ,th6.., improvement  in  beam  quality  attainable  by 
•idaptive  optical  systems  that  introduce  a compensating  phase 
h,i»tiibut;ion  which  .tends  to  cancel  the:  errors  introduced  by 
the*  phase:  aberrations.  The  compensating  phase  distribution 
' }.n  assuinud  to  be  generated  ay  a corrector  *i  rror  located  in  ^ 
tneVVphJ  sm  tr.iim  e.  cj,  , a segwor-ted  dr  Vm'o  unable  mirror 
ItMhtajifl  Just  before  th^  oitpu*.  telesr.-o^ik.  Furthermore',  wc  ■ 
•*r:«  '.yir.f.Spj rued  primarily  v/i<h  close.-  ’,■■* Loop  wstmmn  that  obtain 
1 t.h<>  < Vh-isi-t.,;  jphas a information  via  the  on  I’/rg./  returned  from 
the  t stfi.'ei'.  ( 

. r>he  principal  objective  of  these  studies  has  been  to 
determine  the  effectiveness  of  closed-loop  adaptive  optical 
systems  in  compensating  for  the  effects  of  turbulence  and 
thermal  blooming.  The  effects  of  device  distortion  and 
optical  train  errors  were  not  considered  but  could  easily 
be  included  in  the  formalism  developed  hero  if  the  nature  ot 
these  errors  was  specified  (e.g.,  the  spatial  and  temporal 
spectrum  of  the  errors) . 

To  obtain  quantitative  estimates  of  the  degree  of 
compensation  achievable  with  closed-loop  adaptive  optical 
systems,  we  have  developed  a computer  model  of  such  systems, 
The  model  includes  a servo  system  simulator,  a corrector 
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mirror  si;i)dj  ator , a propagation  simulator  and  a target 
simulator.  The  propagation  simulator  utilizes  a propagation 
algorithm  that  accurately  models  the  optical  distortions 
introduced  by  atmospheric  turbulence  and  thermal  blooming. 

Before  proceeding  to  the  discussion  of  the  simulation 
and  the  results  obtained  therefrom,  we  discuss  in  Section  II 
some  of  the  fundamentals  of  adaptive  phase  compensation  for 
the  effects  of  medium  inhomogeneicies.  This  is  followed  in 
Section  HI  by  a discussion  of  the  detail*  of  the  computer 
simulation  that  has  been  developed  on  this  contract.  Results 
obtained  from  the  computer  simulation  and  from  analysis  are 
th»n  given  in  Section  IV.  The  conclusions  that  we  have 
reached  based  on  this  work  and  our  recommendations  for 
Additional  work  on  the  computer  simulation  of  closed-loop 
adaptive  optical  systems  are  discussed  ?’n  Section  V. 
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11.  / FUNDAMENTALS  OF  ADAPTIVE  COMPENSATION  FOR 

MEDIUM  lNHOMOGENEl'HES 

/ 

/ 

/ In  this  section  vie  discuss  the  basic  principles  ar.d 

/ 

■'  fundamental  limitations  jf  two  Inroad  classes  of  adaptive  optical 
techniques  tor  compensating  v-ai  medium  inhomogenei ties  caused  by 
thermal  blooming  and  turbulence.  The  two  classes  of  techniques 
to  be  considered  are  (ID  return-wave  algorithms  in  which  the 
phase  corrections  a ra  determined  from  direct  measurements  of 
the  returned  field  (e.g.,  this  phase  conjugate  technique)  or 
by  a maximisation  of  some  functional  of  the  irradiarice 
distribution  obtained  by  imaging  the  target  (e.g.,  maximization 
of  the  energy  through  a pinhole  for  a point  glint  target  or 
maximizing  the  ’sharpness"  of  an  image  of  an  extended  target) 
and  (2)  outgoing-wave  multidither  algorithms  In  which  the  phase 
corrections  arc  determined  by  maximizing  the  integrated 
irrndianee  of  the  field  reflected  from  the  target  and  inter- 
cepted by  a tin  tec.  tor  located  in  th  ■?  vicinity  of  the  t ran  a.  mi  tt  ur 
(o.u,  . the  s inuin'ii  i.i- ill  multidither  technique  used  in  .)  0- 
vieme  it  oxo»!r  inv.Mt.a  l COAT  ;.ys‘ em  renstruottn  by  tfTtti) . :cti 
compar  tn<i  those  tec  vnic- .ies  ’,e  w:J.  •.  assume  that  v.h<»  ?»c .a,.-t  ivu 
system  responds  instantaneously  ts  changes  in  1: no  stiite  of 
the  medium  caused  by  thermal  blooming  and  turbulence,  in 
effect,  this  assumption  linearizes  the  propagation  phenomenon 
since  we  then  do  not  have  to  be  concerned  with  the  local 
effect  (in  time)  of  the  phase  correction  on  the  medium.  This 
idealization  should  apply  as  long  as  the  response  time  of  the 
adaptive  system  is  short  compared  to  the  time  constants  for 
the  thermal  blooming. 

For  each  class  of  technique  we  are  interested  in  deter- 
mining the  cffoct  of  the  phase  corrections  on  the  field  at 
the  target.  This  field  cart  be  derived,  in  general,  from  a 
scalar  function  u which  satisfies  the  relation 

u(x~)  - / dx/  u0(x‘)  fl<x~,x‘>  (2.1) 
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where  ug  in  the  field  at  the  transmitter  aperture  and  G is 
the  Green's  function  for  the  inhomogeneous  medium;  i.e., 

G ( k" , x 1 ) is  a scalar  function  from  which  we  can  derive  the 
field  at  x4*  excited  by  a point  source  located  at  x*.  It 
follows  from  Eq.  (2.1)  that  the  field  at  a point  x caused  by  the 
reflection  of  the  transmitted  wave  from  a point  target  at  yr 
is 


Vs> 


G (x > >0  J dx'.  ug  (x* 


) G(jT,X') 


(2.2) 


Likewise,  the  field  at  x caused  by  the  reflection  of  the 
transmitted  wave  from  an  extended  target, having  the  complex 
reflectivity  function  p(x~)  is 

uR(x)  » / dx*“  p(x")  G (x  fX"  > (dx'  (x " ) G(>T  jX')  (2. 3) 


Equations  (2 .2)  and  (2.3)  are  the  fundamental  relations  upon  which 
w-a  sht-11  base  the  discussion  of  the  two  c.larses  a!  adaptive 


camper)  nation  techniques , 
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;v)i  retufn-WAV'f . ihXgcrit'liiia  have  ins.  property  that  thi 
phase  correction  applied  to  the  t.ransm.ictnd  wave  is  the 
conjugate  o Jf  the  phuse  of  t.he  field  returned  from  the  target. 
Thin  if  true  not  only  for  the  explicit  phase  conjugate  technique 
but  also  for  the  case  in  which  a functional  of  the  irradiance 
distribution  obtained  by  imaging  the  target  is  maximised. 

Muller  and  Buffington  have  proved  this  for  an  incoherent  imaging 
system  that  maximizes  the  "sharpness"  of  an  image  of  an  extended 
target.1  It  is  also  easily  proved  for  a coherent  system  that 
maximizes  the  energy  from  a point  glint  target  through  a 
pinhole  in  the  image  plane  of  a receiving  aperture.  In  this 
case  the  field  at  the  receiving  aperture  ie  proportional  to 
G (x  ,x~ ) which  can  be  expressed  in  the  form 
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G(x,x")  « A(x,x~)  exp 


|i <()  (x,x/)  j 


G0(x,V-) 


(2.4) 


where  GQ  ie  the  free  space  Green's  function  and  A and  <{>  are 
the  amplitude  and  phase  variations  caused  by  the  medium 
inhomogeneities.  After  the  introduction  of  a compensating 
phase  distribution  ^ c ( x ) in  the  aperture  plane,  the  field  in 
the  image  plane  of  a lens  used  to  image  the  point  x~  is 


“image <2t> 


/ 


dx  A(x,X“)  exp 


[i (x,x~ 


) + 


(2.5) 


where  f is  the  distance  from  the  aperture  plane  to  the  image 
plane  of  the  len3.  Hence  the ’ irradiance  through  a pinhole 
located  at  x^  * 0 is 


6) 

and  we  note  that  since  A is  positive  the  maximum  value  of 
^ima  is  attained  under  ^ , i.e.,  when  the  compensating 

phaau  distribution  is  the  conjugate  of  the  atmospheric  phase 
distortions . 

Another  important  general  property  of  return-wave 
algorithms  iB  that  they  are  based  on  the  fact  that  Green's 
function  G(x,x")  satisfies  the  reciprocity  relation  (see 
Appendix  A for  a discussion  of  reciprocity) 


1 image  ^ 


I dx,  A(x,x')  exp  £i  ($  <x,x")  + <j)J(x)jj 


(2, 


G (x / x " ) = fl(x',x) 


(2.7) 


which  moans  that  the  field  observed  at  a point  x from  a source 
at  a point  x"  is  the  same  as  the*  field  observed  at  x‘  from 
a source  at  x.  Equation  (2.7)  applies  even  though  the  intervening 
space  between  x and  x'  is  inhomogeneous.  The  only  restriction 
is  that  we  assume  that  the  propagation  is  linear,  which  it 
is  under  the  assumption  we  have  made. 


Let  us  now  consider  the  effect  of  introducing  as  a 
phase  correction  the  conjugate  of  the  phase  of  the  field 
returned  from  the  target.  For  a point  target  we  note  from 
Eq.  (2.2)  that  the  returned  field  ia  proportional  to  G(x,x") 
and  thus  that  the  phase  correction  applied  by  a return-wave 
system  is  The  field  at  the  target  is  then  given  by 


u(x")  - 


/dx'  u9(x')  exp i 0 G(xT,x') 

~ J **'  u8(x')  A(x",x')  G0(x",x')  , 


(2.8) 


where  we  have  used  Eq.  (2.4)  to  express  G in  terms  of  A,<p  and 
the  free  apace  Green's  function  GQ  and  have  used  the  reciprocity 
relation  given  inEq.  (2.7)  to  sliminate  the  phase  terms  <Mx'»>0 
and  <P  (x"  ,x ')  , i . a.  , 


4>(x',#x")  + <t>(x" ,x')  * 0. 


(2.9) 


The  result  inEq.  (2.8)  indicates  that  return-wave  algorithms 
yield  good  compensation  in  the  case  of  point  glint  targets 
as  long  as  the  amplitude  factor  A does  not  vary  significantly 
across  the  transmitting  aperture.  Appreciable  variation  in  A 
across  the  aperture  will  lead  to  an  effective  apodization  of 
the  aperture  and  will  tend  to  increase  the  sidelobe  level 
and  the  beamwidth  of  the  transmitted  energy. 

For  an  extended  target  the  return  field  is  given  in 
(2.3),  which  we  repeat  here  for  convenience 


Eq 


Up (x)  * j dx"  P(x")  G(x,x")  J dx/  u8  fx')  G(x",x'). 


(2.  3) 

If  the  target  reflectivity  function  is  simply  a collection  of 
zero  width  glints  located  over  an  extended  region,  it  can  be 
shown  that  the  previous  point  target  analysis  applies  and  we 
conclude  that  a return-wave  system  will  yield  good  compensation 
for  such  targets  as  long  as  the  amplitude  factor  A does  not 
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vary  too  much.  More  realistically,  however,  the  target  glints 
have  finite  extent  and  the  previous  analysis  does  not  apply. 

For  example,  consider  the  case  where  there  is  a single  glint  of 
finite  extent.  Moreover,  let  us  assume  that  the  range  of  p(xf~) 
is  small  compared  to  that  of  the  incident  field  from  the  trans- 
mitter. In  this  case  the  return  field  can  be  approximated  by 


uR(x) 


p(x~)  G(X 


>X ‘)j  dx'  uB(x')  G(x£,x')\  (2. 


10) 


where  x£  is  the  point  where  the  reflectivity  function  r>(/') 
has  its  maximum  value.  Hence,  a return-wave  algorithm  will 
apply  as  a correction  the  conjugate  of  the  phase  of  the  factor 


i>  (x")  G ( x , x" ) 


(2.11) 


which  is  a weighted  average  of  the  Green's  function  over  the 
domain  of  the  reflectivity  function  c(x").  This  phase  correction 
policy  will  yield  good  compensation  only  if  the  phase  of  G(x,x") 
does  not  vary  appreciably  over  the  range  of  p.  Otherwise,  the 
conjugate  phase  is  garbled  by  the  interference  between  the 
returns  from  the  different  parts  of  the  domain  of  u ( x" ) . The 
negative  of  the  phase  of  the  factor  in  Eq.  (2.11)  will  then 
not  exactly  cancel  the  phase  <t>  (x^f>0  of  G(x“,x‘)  and  the 
compensation  will  be  degraded.  An  example  of  this  type  of 
problem  is  provided  by  a similar  situation  in  adaptive  imaging 
through  inhomogeneous  media.  It  is  well  known  that  good 
compensated  imaging  can  be  obtained  only  if  the  object  lies 
within  an  iBoplanatic  region,  which  has  the  property  that  the 
phase  distortion  <i*  ( x , x“ ) introduced  on  a point  source  field 
radiated  from  a point  x"  on  the  object  to  a point  x in  the 
receiving  aperture  is  nearly  the  same  for  points  x"  within 
this  region. 

The  importance  of  the  Isoplanatic  problem  in  return- 
wave  adaptive  compensation  systems  will  depend  largely  on  the 
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nature  of  the  target  reflectivity  function  and  on  the  source 
of  the  medium  inhomogeneity.  In  our  computer  simulations  of 

phase  conjugate  compensation  for  the  effects  of  thermal 
blooming  we  have  assumed  that  the  target  glints  extend  over  an 
area  approximately  equal  to  the  diffraction-limited  size  of 
the  transmitted  beam.  Under  these  conditions/  we  have  found 
that  the  compensation  achieved  tends  to  be  disappointingly 
small  and  we  attribute  this  to  the  above  discussed  isoplanatic 
problem.  Numerical  problems  associated  with  the  attendant 
wide-angle  scattering  have  prevented  us  thus  far  from  simulating 
glints  that  are  significantly  smaller  than  the  diffraction- 
limited  spot. 

For  turbulence  compensation  problems  we  can  readily 
estimate  the  size  of  the  isoplanatic  region.  Its  linear  extent 
Diso  rou9hlV  ecJual-  tn  the  distance  between  two  points  on 
the  target  at  which  the  mean  square  difference  in  the  phase  of 
a point  source  field  radiated  from  a point  x in  the  transmitter 
aperture  is  unity.  Assuming  Kolmogorov  turbulence/  we  obtain 
the  relation 

Di.0  - (1'1  ci5  f>'3/5  <2'l2> 

2 

where  Cf,  is  tho  index  structure  constant/  k is  the  wavenumber 

2fi / a and  f is  the  distance  to  the  target.  In  the  computer 

simulations  of  phase  conjugate  compensation  for  the  effects 

of  turbulence,  we  have  generally  obtained  good  correction 

with  glints  that  extend  over  a diffraction-limited  spot  area. 

This  is  consistent  with  the  fact  that  the  isoplanatic  region 

associated  with  the  turbulent  Inhomogeneltles  was  larger 

than  the  diffraction-limited  spot  in  these  runs.  For  instance, 

tho  diameter  of  the  diffraction-limited  spot  obtained  with 

a 70  cm,  10.6  urn  laser  beam  focussed  at  a range  of  2 km  is 

approximately  3 cm  (t\/D) , whereas  the  diameter  of  the 

isoplanatic  region  obtained  from  Eq.  (2.12)  is  in  the  range 

6 • D.  ••  26  cm  for  the  range  of  index  structure  constant 
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values  used  in  the  simulations  (5  x 10* 
6 x 10'15  cm"2/3) . 
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2 . 2 Outgoing-Wave  Multidlthor  Algorithms 

Outgoing-wave  multidither  algorithms  adjust  the  phase  of 
the  transmitted  wave  to  maximize  the  irradiance  of  the  field 
reflected  from  the  target.  To  accomplish  this,  the  phase  of 
the  transmitted  field  iB  dithered  in  a distinct  set  of 
spatial-temporal  modes.  The  resulting  modulations  on  the 
irradiance  of  the  field  reflected  from  the  target  are  then 
detected  and  sent  to  a servo  Bystem  that  tends  to  set  the 
corrector  mirror  actuators  so  as  to  drive  the  modulations  to 
zero.  In  doing  so,  the  servo  system  sets  the  actuators  in 
the  positions  where  the  derivatives  of  the  detected  signal 
with  respect  to  the  phase  correction  are  zero,  which  is  the 
condition  for  a maximum  (the  fact  that  it  is  a maximum  and 
not  a minimum  can  be  determined  by  the  phase  of  the  signal). 

The  detector  current  in  an  outgoing-wave  multidither 
system  is  proportional  to  integral  I of  the  irradiance  of  the 
field  uR  reflected  from  the  target.  For  a point  target,  we 
obtain  from  Eq.  (2.2) 


I ■ c / dx  | G(x,x“)  I 2 (2  . 13) 

where  c is  the  irradiance  of  the  transmitted  field  at  the 
target 


c *=  | / dx'  uB(x‘)  G(x,x‘)|2  (2.14) 

Note  that  the  integral  in  Eq.  (2.13)  does  not  depend  on  tne 
transmitter  field  u . Hence,  a maximization  of  I with  respect 
to  the  transmitter  phaBe  yields  maximum  irradiance  at  the 
target,  as  desired. 


For  extended  targets  the  reflected  field  is  given  by 
Eq.  (2.13).  as  in  the  case  of  return-wave  algorithms,  the  above 

point  target  analysis  applies  if  the  target  reflectivity 
function  is  simply  a collection  of  zero  width  glints  located 
over  an  extended  region  (with  the  difference  that  the  receiver 
detector  must  now  also  be  extended  in  order  to  achieve  the 
desired  maximization  of  the  irradiance  of  the  target).  An 
an  example  of  an  extended  target,  we  shall  consider  the  case 
discussed  in  the  return-wave  algorithm  discussion  of  Section 
2.1;  viz.,  a single  glint  of  finite  extent  for  which  the  range 
of  the  reflectivity  function  p(x")  is  small  compared  to  that 
of  the  incident  field  from  the  transmitter.  In  this  c*ae  the 
return  field  can  be  expressed  in  the  form  given  in  Eq.  (10) 

A, 

and  the  corresponding  expression  for  tho  integrated  irradiance 
returned  from  the  target  becomes 

I - c / dx  |/  dx"  p(x“)  G(x,x")|2  (2.15) 

where  c 1b  again  the  irradiance  of  the  transmitted  field  at  the 
target.  Although  the  field  returned  from  the  extended  glint 
differs  from  that  returned  from  the  point  glint,  we  note  from 
Eq.(2.15)  that  a maximization  of  I again  leads  to  a maximization 
of  the  irradiance  at  tho  target.  As  in  the  point  target  case, 
this  occurs  because  tho  integral  over  x in  Eq.  (2.15)  does  not 
depend  on  the  transmitter  field  u . 

Tho  above  analysis  indicates  that,  in  contrast  to 
return-wave  algorithms,  outgoing-wave  multidithor  algorithms 
are  relatively  insensitive  to  target  glint  size.  Of  course, 
this  is  true  only  up  to  a point.  If  the  range  of  the  glint 
reflectivity  function  p(x")  is  large  compared  to  that  of  the 
field  incident  from  tho  transmitter,  the  return  field  is 
approximately  given  by 

uR(x)  - p(>£)  / dx"  G(x,x")  / dx*  ug(x*)  G(x",x*),  (2.16) 
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and  the  integrated  irradiance  returned  from  the  target 
becomes 

1 “ |p<*£)  I2  / **  1/  d2f  0(x,>T ) ; dx  uM(x')  C3  (^'  ) I 2 

(2.17) 

Hence,  in  this  case  an  outgoing-wave  multidither  ayntem 
maximizes  a complicated  functional  of  the  transmitter  irradi- 
ance at  the  target  instead  of  the  transmitter  irradiance 
itself  and  we  should  not  expect  to  achieve  good  compensation 
(this  behavior  has  been  observed  experimentally  with  the 
HRL  10-element  multidither  COAT  system). 

In  summary,  the  results  of  this  section  and  those  of 
Section  2.1  indicate  that!  (1)  both  return-wave  and  outgoing- 
wave  (multidither)  phase  compensation  algorithms  should  yield 
good  compensation  for  medium  inhomogeneities  for  point  glint 
targets;  (2)  return-wavs  algorithms  are  much  more  sensitive 
to  glint  size  than  outgoing-wave  algorithms,  especially  in 
the  cane  of  thermal  blooming  compensation;  and  (3)  when  the 
size  of  the  target  glints  significantly  exceeds  the  size  of 
the  diffraction-limited  focal  spot  of  the  transmitted  field 
at  the  target,  we  do  not  expect  the  outgoing-wave  multidither 
algorithm  to  work. 
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Til. 


ADAPTIVE  OPTICAL  SYSTEM  SIMULATION 


3.1  Organisation  and  Approach 


The  essential  elements  of  the  adaptive  optics  computer 
simulation  are  indicated  in  Fig.  3.1.  The  simulation  contains 
computer  routines  that  model  the  operation  of  an  adaptive 
optical  system  consisting  of  a servo  system  and  a corrector 
mirror.  The  effect  of  the  atmosphere  on  the  resulting  laser 
beam  propagation  is  modeled  by  a tima-dependent  propagation 
code  that  treats  the  effects  of  both  turbulence  and  thermal 
blooming.  The  simulation  models  a closed-loop  adaptive  system 
in  which  the  control  information  is  provided  by  the  detection 
of  a return-wave  reflected  from  the  target.  Propagation  on  the 
return  path  is  modeled  by  a time-dependent,  linear  propagation 
code,  which  includes  the  effects  of  tho  inhomogeneities 
caused  by  turbulence  and  the  absorptive  heating  introduced 
by  the  high-power  transmitted  wave.  The  nature  and  current 
status  of  the  various  parts  of  the  simulat ion  are  discussed 
in  Sections  1.2  to  3.5. 


i 

t 
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As  indicated  above,  we  use  a ti me-dopcndcnt  propagation 
code  to  model  tho  laser  propagation  on  the  outgoing  and  return 
paths.  This  differs  from  tho  approach  taken  by  Herrmann  and 
Bradley  of  Lincoln  Labor afar  y who  h.ivr  1. 1. c d u c«*  .. .... . * 

their  phase  compensation  studies.  It.  is  admittedly  true  that 
the  use  of  a time-dependent,  propagat  i on  code  greatly  increases 
the  complexity  and  cost  of  the  calculation  but:  we  have  felt 
compelled  to  use  this  approach  for  two  reasons.  First,  the 
time-dependent  approach  it  closer  to  the  physics  of  the  stri- 
atton  that  we  are  trying  to  simulate  since  we  are  primarily 
interested  in  adaptive  optical  systems  that,  respond  in  a time 
short  compared  to  the  Lime  it  takes  the  atmosphere  to  change. 

Huch  adaptive  systems  perform  an  essentially  instantaneous  deter- 
mination of  tho  required  phase  front  . Moreover,  the  dotormi nation 
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Fig.  3.1.  Adaptive  optical  system  simulation  flow  diagram. 
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is  continuously  updated  as  the  medium  changes,  rather  that  to 

be  put  in  and  held  fixed  until  a steady  state  is  reached,  as 
is  implicitly  assumed  whan  a cw  propagation  code  is  used. 

The  distinction  between  the  time-dependent  and  cw  approaches 
is  especially  important  for  the  simulation  of  adaptive  systems 
that  employ  return-wave  algorithms.  Whereas  it  can  be  argued 
that  a cw  code  in  conjunction  with  a maximization  routine 
closely  models  the  behavior  of  an  adaptive  optical  system 
employing  an  outgoing-wave  multidither  algorithm,  it  does  not 
appear  possible  to  model  return-wave  systems  with  the  cw  approach 
since  such  systems  do  not  explicitly  maximize  the  irradiance  on 
the  target,  instead,  as  discussed  in  Section  2.1,  they  continu- 
ously sec  the  phase  of  the  transmitted  wave  equal  to  the  conjugate 
of  the  phase  of  the  field  returned  from  the  target.  There  does 
not  appear  to  be  any  reason  to  expect  that  the  final  state 
reached  by  such  a system  can  be  obtained  from  a cw  calculation 
or  a sequence  of  such  calculations. 

The  other  reason  for  the  use  of  a time-dependent  propa- 
gation code  in  our  studies  is  the  desire  to  follow  the  time 
development  of  the  correction  process,  especially  in  the  case 
of  adaptive  systems  employing  outgoing-wave  multidither  control 
algorithms.  Ab  discussed  in  Section  2.2,  an  outgoing-wave 
multidither  control  system  attempts  to  maximize  the  irradiance 
at  the  target  at  each  instant.  We  wish  to  know  if  this 
instantaneous  maximum  control  policy  leads  to  a global  maximum 
in  the  thermal  blooming  correction  problem  or,  instead,  yields 
a secondary  maximum.  This  can  be  determined  by  comparing  the 
maximum  obtained  with  the  time-dependent  and  the  cw  codes 
because  the  cw  approach  yields  the  global  maximum. 

The  time-dependent  propagation  algorithm  used  in  the 
simulation  is  basically  a multi-pulse  code  in  which  the  laser 
energy  is  propagated  in  a series  of  short  pulses.  The  pulse 
length  is  assumed  short  enough  that  no  heating  occurs  during 
the  pulse.  Furthermore,  it  is  assumed  that  the  heating 
caused  by  a pulse  has  reached  a steady  state  by  the  time  the 
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next  pulse  arrives.  The  term  i.teady  state  is  used  here  in  the 
sense  that  all  density  gradients  traveling  at  the  sound  velocity 
are  assumed  to  have,  left  the  region  occupied  by  the  laser  beam 
by  the  time  the  next  pulse  arrives.  The  only  gradients  present 
are  those  that  are  convected  with  the  local  crosswind  velocity. 
Each  pulse  "sees"  the  medium  created  by  the  steady  state  heating 
of  all  previous  pulses  but  is  unaffected  by  its  own  heating. 

This  model  allows  us  to  simulate  either  pulsed  or  cv;  propagation, 
depending  on  the  time  between  pulses  relative  to  the  time  it 
takes  an  air  parcel  to  traverse  the  beam.  CW  propagation 
conditions  are  simulated  by  allowing  significant  overlap  between 
the  regions  heated  by  successive  pulses. 

The  time  sequence  of  events  in  the  simulation  is  shown 
An  Fig.  3.2.  First,  a high-power  pulse  of  laser  energy  is 
propagated  to  the  target  and  the  resulting  change  in  the 
density  of  the  air  in  the  region  occupied  by  the  beam  is 
calculated,  taking  into  account  the  convective  effect  of  the 
local  crosswind  velocity.  It  is  assumed  that  the  adaptive 
phase  compensation  is  applied  just  prior  to  the  firing  of  the 
next,  high-power  laser  pulse.  Moreover,  we  are  presently 
assuming  that  the  adaptive  optical  system  has  infinite  band- 
width bo  that  it  is  able  to  determine  and  apply  the  phase 
correction  instantaneously  without  any  delay  or  time-averaging 
effects.  The  information  required  to  determine  the  phase 
correction  is  obtained  by  transmitting  a low-power  reference 
pulse  to  the  target  and  reflecting  the  same  from  the  target. 

This  reference  pulse  propagates  through  tha  medium  that  will 
be  seen  by  the  next  high-power  pulse  and  thus  provides  the 
necessary  phase  correction  information.  This  process  is  then 
repeated  for  each  subsequent  nigh-power  pulse. 

With  the  exception  of  some  runs  in  which  the  effect 
of  the  atmosphere  was  represented  by  a Bingle  phase  screen 
in  front  of  the  transmitter,  all  of  the  data  that  have  obtained 
to  this  date  apply  to  the  simulation  of  adaptive  systems 
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utilizing  a return-wave  control  algorithm,  i.e.,  systems  that 
correct  by  inserting  the  conjugate  of  the  phase  of  the  wave 

returned  from  the  target.  The  reason  for  this  concentration 
on  the  phase  conjugate  simulation  is  simply  that  such  simu- 
lations are  significantly  easier  to  implement  and  less  costly 
to  run  than  outgoing-wave  simulations.  The  simulation  of 
a return-wave  system  requires  simply  that  we  transmit  a 
single  low-power  ref erence  wave  to  the  target  and  a single 
low-power  reflected  wave  back  to  the  transmitter  aperture. 

The  phase  correction  is  then  equated  to  the  conjugate  of  the 
phase  of  the  reflected  field.  In  contrast,  an  accurate 
multidither  simulation  requires  the  sending  and  receiving  of 
many  low-power  reference  pulses  for  each  high-power  pulse. 

We  must  sample  the  modulations  of  the  integrated  irradiance 
of  the  target  returns  at  a rate  set  by  the  highest  dither 
frequency.  For  example,  it  was  found  in  our  earlier  multi- 
dither servo  system  simulations  that  a sampling  time  increment 
of  5 i.j s e c is  required  for  an  18-channel  system  having  an 
upper  dither  frequency  of  32  kHz  and  an  overall  convergence 
time  of  1 ms.3  Hence,  approximately  200  samples/convor jence 
time  were  required  in  those  simulations.  It  was  practical 
to  take  200  samples/convergence  time  in  the  earlier  work 
because  of  the  simplicity  of  the  propagation  model  that  was 
used.  In  the  present  case,  however,  the  propagation  model 
is  quite  complicated  and  relatively  costly  to  run  and  thus 
it  is  clearly  not  practical  to  do  an  equivalent  simulation. 

At  the  outset  of  this  contract  we  had  planned  to 
circumvent  the  multidither  simulation  problem  as  follows j 

(1)  propagate  a low-power  reference  pulse  to  the  target; 

(2)  propagate  the  associated  low-power  pulse  reflected  from 
the  target  back  to  the  transmitter  aperture;  (3)  use  the 
calculated  return  wave  to  estimate  the  Green's  function  of 
the  medium;  (4)  dither  the  phase  of  the  reference  pulse  and 
calculate  the  resulting  variations  in  the  irradiance  received 


18 


from  the  target,  utilizing  the  Green's  function  information 
obtained  previously  (i.e.,  we  did  not  plan  to  use  the  full 
propagation  simulation  every  dither  cycle,  but  rather  intended 
to  do  this  only  once  per  high-power  pulse);  and  (5)  process 
the  variations  in  received  irradiance  in  a multidither  servo 
simulator  to  determine  the  mirror  actuator  settings.  During 
the  course  of  the  work  on  this  contract,  it  became  apparent, 
however,  that  such  a simulation  1b  an  empty  exercise  because 
it  does  not  simulate  the  desired  behavior  of  an  actual  outgoing- 
wave  multidither  system  in  a thermal  blooming  situation.  The 
problem  is  that  the  assumption  that  the  Green's  function  can 
be  determined  from  the  return  field  reflected  from  the  target 
is  not  valid  in  the  thermal  blooming  problem  unless  the 
target  is  extremely  small,  in  fact,  unless  it  is  nearly  a 
point  target.  But  it  can  be  shown  from  the  analysis  in 
Section  II  that  outgoing-wave  multidither  and  return-wave 
adaptive  optical  systems  both  put  in  the  conjugate  phase  when 
the  target  is  a point  target.  Therefore,  the  simplified 
approach  to  multidither  outlined  above  either  yields  results 
equivalent  to  those  obtained  with  a return-v,ave  algorithm 
(point  glint)  or  yields  results  that  are  patently  incorrect 
(extended  glint) . 

We  now  think  that  the  beat  way  to  simulate  outgoing 
multidither  adaptive  systems  is  to  use  a function  maximization 
routine  in  conjunction  with  the  time-dependent  propagation 
code.  This  is  a modified  version  of  the  approach  used  by 
Herrmann  and  Bradley  in  their  phase  compensation  studies. 

They  have  used  a function  maximization  routine  in  conjunction 
with  a cw  propagation  code.  We  recommend  the  use  of  a time- 
dependent  code  because  it  will  allow  us  to  determine  directly 
whether  or  not  the  outgoing-wave  multidither  algorithm  yields 
a truly  global  maximum  and  not  juBt  a secondary  maximum. 
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As  discussed  in  Section  3.1,  we  have  primarily  used  a 
phase  conjugate  control  policy  in  the  computer  studios  of 
thermal  blooming  and  turbulence  compensation  performed  on 
this  contract.  Moreover,  the  phase  conjugate  control  algorithm 
that  we  have  used  is  particularly  simple  in  that  wa  have 
simply  equated  the  phase  of  the  transmitted  wave  to  the  conjugate 
of  the  phase  of  the  wave  returned  from  the  target  on  a mesh 
point  by  mesh  point  basis.  We  have  not  taken  into  account 
the  effect  of  the  finite  area  of  the  phase  detector  and 
detector  noise  in  these  studies  because  we  wiBhed  to  establish 
an  upper  limit  on  the  amount  of  correction  that  is  achievable 
with  a phase  conjugate  control  algorithm  and  to  avoid  the 
problem  of  determining  whether  the  particular  detector  configu- 
ration used  is  the  optimum  choice.  For  example,  there  is  the 
question  of  how  large  to  make  the  detectors.  Should  we  use  a 
collection  of  small  point  detcsctors  or  should  wa  ubs  larger 
detectors  that  yield  spatially  averaged  phase  estimators?  The 
refinements  in  our  present  approach  necessary  to  account  for 
these  effects  are  discussed  below. 

The  phase  detector  can  be  represented  as  a finite 
aperture  heterodyne  detector,  which  performs  a coherent 
summation  of  the  return-wave  field  over  the  detector  aperture, 
followed  by  a device  that  determines  the  phase  of  the  detector 
output.  The  detector  output  is 

idet  a J d*  UR(*>  uIo(S>  + inoise  t3<1) 

aperture 

where  p is  a constant,  uR  is  the  field  returned  from  the 
target,  uL0  is  the  local  oscillator  field,  and  inoise  Is  a 
shot  noise  current.  Note  that  even  in  the  absence  of  noise, 
a finite  detector  does  not  yield  an  exact  measure  of  the 


phase  of  the  return-wave  field  but  rather  gives  an  average 
value  determined  by  the  integral  in  Eg.  (3.1).  The  phase  of 
the  detector  output  is 


where 


pdet 


■ tan 


-1 


An  sind 

tan,V  + jq 

.AH" 

As  5oi^ 


N 

S 


r. 


(3.2) 


(3.3) 


^ noise  - AN  exP(i  V 


For  shot  noise-limited  detection,  the  noise  current  in  Eq. 

(3.4)  is  a gaussian  random  process  and  thus  can  be  simulated 
by  a gaussian  random  number  generator.  The  mean  square 
magnitude  of  these  numbers  is  determined  by  the  detector 
signal-to-noise  ratio  S/N 


- Ag  (S/N)  *”1  (3.5) 

For  each  return  wave  a finite  set  of  phase  samples 
can  be  determined  from  the  relations  given  in  Eqs.  (3.2)  through 

(3.5).  The  conjugate  of  three  samples  can  then  be  used  as 
inputs  to  the  mirror  simulator  discussed  in  Section  3.3.  The 
resulting  phase-corrected  wave  is  then  transmitted  to  the 
target  and  a new  return  wave  calculated.  For  point  targets, 
the  necessary  phase  correction  is  established  by  a single 
calculation  of  this  type.  Our  experience  indicates  that  this 
transmit-return  wave  calculation  must  be  iterated  four  to  five 
times  in  order  to  achieve  convergence  on  a complex  target, 
i.e.  , a target  having  more  than  one  glint. 
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With  regard  to  the  question  of  how  to  simulate  outgoing- 

wave  multidither  control  algorithms  we  recommend,  as  mentioned 
in  Section  3.1,  that  a function  maximization  routine  be  used 
in  conjunction  with  a time-dependent  propagation  code  (see 
Fig.  3.3  . The  basis  of  this  approach  is  the  fact  that  an 
outgoing-wave  multidither  algorithm  is  a maximization  process. 

It  attempts  to  maximize  the  integrated  irradiance  of  the  field 
returned  from  the  target.  Hence  the  detector  in  Fig,  3.3  is 
represented  in  the  simulation  by  an  integration  of  the  return 
irradiance  over  an  area  equal  to  that  of  the  receiver  being 
modeled.  The  maximization  can  be  done  with  respect  to  any 
desired  set  of  parameters  that  describe  the  state  of  the 
corrector  mirror,  e.g.,  the  actuator  position,  the  coefficients 
of  a set  of  orthogonal  polynomial  such  as  Zernike  functions,  or 
the  coefficients  of  some  predicted  set  of  correction  functions 
(such  as  parabolic  refocuB  adjustment  plus  a variable  percentage 
of  the  phase  distortion  introduced  in  the  first  propagation 
step  as  implemented  by  Bradley  and  Herrmann) . 

3,3  Mirror  Simulation 


We  have  developed  simulations  of  both  segmented  and 
dei'ormable  mirrors.  The  segmented  mirror  simulation  is  a 
relatively  trivial  exercise  in  that  it  simply  consists  in 
dividing  the  mirror  mesh  into  a set  of  square-shaped  zones 
within  an  outer  circular  boundary  (see  Fig.  3.4).  The  phase 
within  each  zone  is  assumed  to  be  constant;  i.e,,  the  mirror 
consists  of  a set  of  flat  segments. 

In  the  deformable  mirror  simulation,  it  is  assumed 
that  the  mirror  is  driven  by  actuators  equi-spaced  on  a square 
array  of  mesh  points.  The  fundamental  problem  iB  to  find  a 
relatively  simple,  low-cost  method  for  calculating  the  mirror 
profile  resulting  from  an  arbitrary  set  of  actuator  displace- 
ments. It  is  not  practical  to  use  a sophisticated  simulation 
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such  as  the  NASTRAN  code  in  the3e  studies.  Instead,  we  have 
simplified  the  problem  by  assuming  that  thin  piste  theory  is 
applicable,  and  that  linear  superposition  holds.  It  is  alBO 
assumed  that  the  connection  to  the  mirror  from  each  actuator 
is  effected  at  a single  point  (i.e.,  the  actuators  have  zero 
width  at  the  point  of  contact  with  the  mirror)  and  that  purs 
thrust  without  torque  is  transmitted.  With  these  assumptions, 
the  problem  reduces  to  finding  the  mirror  profile  resulting 
when  one  actuator  is  raised  unit  height  with  all  of  the 
surrounding  actuators  held  fixed  at  zero  height.  This  is 
referred  to  as  the  unit  profile.  The  total  mirror  profile  at 
any  instant  is  calculated  by  summing  a set  of  unit  profiles, 
each  centered  at  one  of  the  actuators  and  each  weighted  by 
the  height  of  its  associated  actuator. 

The  detailed  form  of  the  unit  profile  depends  on  the 
actuator  configuration.  We  have  calculated  unit  profiles  for 
deformable  mirrors  having  37,  57,  and  97  actuators.  In  each 
case  the  unit  profile  is  calculated  as  follows.  First,  we 
consider  the  mirror  profile  that  results  from  the  displacement 
of  a single  actuator  at  the  center  of  a thin  plate  which  is 
constrained  only  at  a circular  edge  whose  radius  is  large 
compared  with  the  extent  of  the  region  occupied  by  the  mirror. 
No  other  actuator  constraints  are  assumed  at  this  stage.  We 
determine  the  profile  from  the  following  relation  obtained 

A 

from  thin  plate  theory 

displacement  « 5(a2  - r2)  - 4r2  in  ~ (3.6) 

where  a is  the  distance  to  the  edge  of  the  plate  and  r is 
the  distance  from  the  central  actuator  to  the  point  at  which 
the  displacement  is  measured. 

For  the  profile  in  Eq.  (3.6)  the  height  at  adjacent 
actuator  positions  is  not  zero  when  the  central  actuator  is 
displaced.  To  determine  the  unit  profile,  we  sum  the 
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contributions  to  the  total  height  at  a given  actuator  position 
from  that  actuator  and  all  surrounding  actuators.  The  height 
of  each  of  the  unconstrained  actuator  displacements  is  taken 
as  an  unknown  coefficient  in  a set  of  simultaneous  equations 
which  are  solved  on  the  basis  that  the  total  height  should  be 
unity  at  the  central  actuator  and  zero  at  the  surrounding 
actuators.  Having  determined  the  unconstrained  actuator 
heights  required  to  give  the  desired  constraint  that  the  net 
displacement  is  zero  at  each  of  the  actuators  except  for  the 
central  one#  where  it  is  unity,  wo  use  Eq.  (3.6)  to  determine 
the  unit  profile  at  all  intervening  mesh  points. 

If  the  profile  is  used  in  the  form  determined  above, 
no  additional  adjustments  need  be  made.  For  convenience, 
however,  we  have  truncated  these  profiles  along  square 
boundaries.  Hence,  we  slightly  adjust  the  resulting  profile 
to  ensure  that  it  satisfies  the  condition  that  if  all  actuators 
are  displaced  at  unit  height,  the  mirror  surface  is  also  at 
unit  height  everywhere.  The  unit  profile  obtained  for  a 
57-actuator  mirror  .is  given  in  Fig.  3.5.  The  unit  profiles 
obtained  for  the  37-  and  97-actuators  mirrors  have  the  same 
general  characteristics,  viz.,  a monoto'nic  decrease  to  zero 
at  the  noarost  neiqhbor  actuators  and  neqative  values  in  the 
space  between  thesa  neighboring  actuators  anct  the  next  set 
of  actuators. 

The  actuator  configurations  utilized  in  the  simulations 
are  shown  in  Figs.  3.6,  3.7  and  3.8.  Note  that  extra  actuators 
arc  placed  around  the  periphery  of  the  mirror.  These  actu- 
ators are  used  to  simulate  a mirror  having  an  unconstrained 
edge.  They  tend  to  Bmooth  the  transition  between  the  region 
occupier!  by  the  actuators  and  that  which  is  beyond. 
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6.  Actuator  configuration  used  in 
37-actuator  deformable  mirror 
simulation  (37  actuators,  20 
dummies,  13  mesh  points  per 
actuator. 
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Actuator  configuration  used  i 
3 /-actuator  deformable  mirror 
simulation  (57  actuators,  24 
dummies,  9 mesh  points  per 
actuator ) , 
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Fig.  3. 


Actuator  configuration  UBed  in 
97-actuator  deformable  mirror 
simulation  (97  actuators,  32 
dummies,  5 mesh  points  per 
actuator) . 


3 . 4 Propagat ion  Simula t io n 


3.4.1  Formulation  — The  field  of  a high-power  laser 
beam  propagating  through  a turbulent  atmosphere  can  be 
derived  from  a scalar  function  u(x)  that  satisfies  the 
equation 

3 2 

i 2 0 + v*  u + k2(n2 (x,t)-l)  u (3.7) 

3 2 

where  is  the  transverse  Laplacian  i)2/3x2  + a2/ay2,  k is 
the  wavenumber  2n/\,  and  n(x, t)  is  the  refractive  index  of 
the  medium.  In  writing  Eq,  (3.7),  we  have  assumed  that  the 
nominal  direction  of  propagation  is  along  the  z axis  of  an 
x,y,z  coordinate  system.  At  laser  wavelengths  the  function 
u(x,t)  varies  much  mora  rapidly  transverse  to  the  direction 
of  propagation  than  it  does  along  the  direction  of  propa- 
gation. This  enables  us  to  make  the  paraxial  approximation 
wherein  the  term  on  the  right  side  of  Eq.  (3.7)  is  neglected. 
We  then  have 

i 2k  ~ + V2  u + k2(n2(x,t)-l)  u * 0 (3.8) 

For  given  initial  data,  i.e.,  values  of  the  field  at 
points  in  the  aperture  of  the  laser  transmitter  (or  on  the 
target),  the  propagation  simulator  must  generate  the  corres- 
ponding field  values  at  the  target  (or  in  the  receiving 
aperture)  such  that  Fq.  (3.8)  is  satisfied.  To  do  thiB,  we 
divide  the  medium  into  slabs  defined  by  planes  nn  which  z is 
constant.  In  going  from  one  slab  to  the  next,  we  write  u(x) 
in  the  form 

u(x,t.)  = cxp(r  (x,t))  i. >(x,t)  (3.9) 

where  I (x, t)  is  a phase  function  associated  with  the  medium 
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inhomogenei ties  between  z-Az/2  and  z+Az/2 


r ( x , t ) 


ik 

2' 


/•  Z+Az/2  / o v 
/ dz'fn  (x,y,z',t)-l) 
J z-Az/2  N ' 


with  the  initial  condition 


(3.10) 


The  complex  amplitude  <ii(x)  then  satisfies  the  equation 

|.i  2k  ^ + exp  ( — T ) V * exp ( r ) ) u = 0 (3.11) 


uj(x,y,z,t)  « u(x,y,z,t> 


(3.12) 


Physically,  these  equations  aporoximate  the  propagation  in 
the  inhomogeneous  medium  by  a two-step  process  at  each  z 
increment.  First,  we  propagate  the  field  u(x,t)  at  z-Az/2 
to  z+Az/2,  assuming  that  the  intervening  space  is  homogeneous. 
The  effect  of  the  inhomogeneities  between  z-Az/2  and  z+Az/2 
is  then  accounted  for  by  multiplying  this  solution  by  the 
phase  factor  exp(l'). 

3.4.2  Converging  Coordinate  Transformation  — • To  reduce 
the  size  of  the  mesh  required  to  solve  Eq.  (3.11)  numerically, 
we  use  converyitv]  coordinates  defined  by  the  transformation 


N (z ) 


(3.13) 


(3.14) 

(3.15) 


(3. 16) 
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fci  ^2  - u — ~k  (3.17) 

where  pQ  is  a characteristic  dimension  of  the  laser  bean  at 
the  transmitter  (e.g.,  the  e-folding  radius  of  a gaussian 
beam)/  f is  the  distance  to  the  focus,  and  u is  a constant 
determined  by  the  requirement  that  the  oolution  be  confined 
within  the  boundaries  of  the  mesh  at  the  focal  plane.  In  the 
absence  of  thermal  blooming  and  turbulence,  the  choice  a - 1 
yields  a coordinate  system  that  converges  at  a rate  determined 
by  the  free-space  diffraction  of  a gaussian  beam  havinq  an 
e-folding  radius  pQ.  With  turbulence  and  thermal  blooming 
effects  present,  we  typically  choose  « in  the  range  1.5  to 
3.0,  depending  on  the  amount  of  additional  spreading  induced 
by  these  effects. 

When  written  in  terms  of  the  converging  coordinate 
variables  defined  above,  Eqs.  (3.9)  and  (3.11)  for  the  complex 
amplitude  are  replaced  by  the  relations 


uj(x,y,z,t)  « £(£.,;)  exp(f)  v(£,r,,t)  (3.18) 


<1*0  ■ (<»1/2 
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(3.19) 

(3.20) 

(3.21) 


where  A r,  1b  the  increment  in  f.  in  going  from  z-az/2  to  z+az/2. 
The  initial  condition  for  v is 


v(^,r,,t)  * ui  (x,y,z,t)/fii  f^,r.) 


(3.22) 
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To  solve  Eq.  (3.21)  we  utilize  the  fact  that  for 

sufficiently  small  values  of  A£  (i.e.,  Az)  the  effect  of  the 
exponential  factors  exp(±f)  in  this  equation  is  smAll.  Hence, 
we  solve  the  simpler  equation  obtained  when  these  factors 
are  equated  to  unity 


v ■ 0 


(3.23) 


3.4.3  Method  of  Solution  — We  use  a fast  Fourier 
transform  technique  to  solve  Eq,  (3.23).  The  basis  of  this 
approach  is  the  fact  that  the  solution  of  Eq.  (3.23)  can  be 
expressed  in  the  form  of  a discrete  Fourier  series 


N-l  N-l 

vU,£,t)  * £ £ t/mn(*'t)  exP|i{Pm5l+qnc2,|  (3l24) 

m-o  n“0 


where  the  Fourier  coefficients  I / are  determined  from  the 

mn 

initial  data  and  Eq.  (3.23)  as  follows.  The  initial  values 
of  1/  are  obtained  by  taking  the  discrete  Fourier  transform 
of  the  initial  values  of  over  a mesh  of  points 

q - ( A-  (N/2 ) ] A r,t  j,2  - ( j - (N/2 ) ] Ar,  U,j  “ 0,1 , . . . ,N-1) 


Wi'^ 


(-1) 


m+n 


N-l  N-l 


i i v((e  - 1)  o - 1) 


exp|-  (mi  + nj)| 


(3.25) 


The  F,  dependence  of  1/  is  then  determined  by  substituting 
Eq.  (3.24)  in  Eq.  (3.23),  which  yields 
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from  which  it  follows  that 


j W*'*'  “ axP 

I : 

f 

> 

Finally , it  can  be  shown  that  in  order  for  the  discrete 
Fourier  series  representation  of  v given  in  Eq.  (3.24)  to  be 
real  when  v is  real,  the  coefficients  pm  and  qn  must  have  the 
form 

Pm  - m (m  - if)  <3-28 

%mm(n’i)  <3'29 

Hence,  for  discrate  points  ■ U - N/2)Ac#  Cj  ■ (j  - N/2)A; 
U,  j - 0,1,. ..  ,N-1) 

v^a  " ?) ( i - r)  r*  » ^ 
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(3.30 


where  B - 2 tt 2 A £/  (A; ) 2 . Note  that  v is  simply  (-1)*+^  times 
the  discrete  Fourier  transform  of  (-l)m+n  V (£,t). 


The  effect  of  the  medium  and  the  factor  exp  [ ( -i/2  ( t. . 

2 1 
+ 1 introduced  by  the  coordinate  transformation  is  taken 

into  account  at  each  step  in  the  calculation  by  multiplying 

the  value  of  v obtained  in  the  previous  step  b.-  the  quantity 

exp(I')  defined  in  Eq.  (3.20),  i.e.,  the  initial  value  inserted 

in  Eq.  (3.25)  is  exp(T)  times  the  value  of  v determined  from 

the  previous  steps. 

3.4.4  Effect  of  Absorption  — The  absorptive  contri- 
bution to  f is  calculated  from  the  relations 


iky_  /■  z+Az/2  /.  \ 

- -r- /E.A,/2  “'(4 


(3.31) 
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(3.33) 


where  y is  the  ratio  of  the  specific  heat  of  air  at  constant 
pressure  to  that  at  constant  volume  (y  = 1.4),  <»  is  the  mole- 
cular absorption  coefficient  of  the  air,  aT  is  the  total 
extinction  coefficient,  vfj  is  the  transverse  wind  velocity 
assumed  to  be  in  the  x direction,  P is  the  ambient  air  pressure 
u is  the  slew  rate  (slewing  assumed  also  to  be  in  the  x direc- 
tion) , E . is  the  energy  per  pulse,  and  y is  the  molecular 

U X HG  O 

polarizability  of  air. 

In  writing  Eq.  (3.32)  we  have  assumed  that  the  laser 
energy  is  delivered  in  a serieB  of  pulses  for  which  the  pulse 
length  is  short  enough  that  no  heating  occurs  during  a pulse. 


Moreover,  we  have  assumed  that  the  interval  between  pulses 

x,  u 

is  long  enough  that  the  heating  caused  by  the  m pulse  is 
fully  developed  by  the  time  that  the  m + 1st  pulse  arrives 
(i.e.,  transient  density  changes  traveling  at  the  velocity  of 
sound  are  assumed  to  have  departed  from  the  region  occupied 
by  the  beam) . 

3.4.5  Effect  of  Turbulence  — The  contribution  to  f 
from  turbulence  is  calculated  by  a Fourier  transform  technique 
which  yields  a two-dimensional  distribution  of  gaussian 
random  numbers  having  correlation  properties  consistent  with 
those  associated  with  the  phase  variations  caused  by  Kolmogorov 
turbulence.  In  particulars,  we  represent  rtur^  by  the  Fourier- 
Stieltjes  integral 

^turb(-}  " 1 f dr{* exp(i  - ‘ x)  (3.34) 

The  increments  df  satisfy  the  relations 
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where  (K)  is  the  spatial  spectrum  of  the  phase  variations  and 
K is  the  spatial  wavenumber.  Note  that,  with  this  representa- 
tion, I' t ur jo  satisfies  the  Weiner-Khintchin  theorem 
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as  required  by  any  valid  representation  of  a homogeneous 
random  process. 

We  have  utilized  a modified  von  Karman  spectrum  for  <t>p 


♦ f (K) 


. 207  C2  k2  Az  K“11/3 
N o 


(3.38) 


where  is  the  index  structure  constant,  Az  is  the  z incre- 
ment, and  K is  2i  divided  by  the  outer  scale  of  turbulence 
Lq.  The  Fourier  integral  in  Eq.  (3.34)  is  evaluated  by  the 
following  discrete  Fourier  transform  approximation  of  this 
integral 


1 turb 
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(3.39) 


where  the  a , ^ are  complex  gaussian  random  numbers  (we  assume 
that  i'  fa  is  gaussian)  with  zero  mean  and  the  following 
correction  properties 


aff.j  = aN-2.,N-j 


(3.40) 
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Real^a^)  Real (u ^ ,)J  - 0 

and  it '(‘N-ii.j  VN-j 

(3.42) 

|lmag(apj)  ImagCa^ ,)])  - 0 

and  tVN-t.J^K-J 

(3.43) 

^Real(a^)  Real (a^ , j ,)^  ■ 1 

or  r-N-Jt.J'-N-l 

(3. 44) 

^Imag(a^ ) ImagCa^  • 1 

a'-fc.j'-J 

(3.45) 

^lmag(a^j)  Imag^ - -1 

» J ‘"N-J 

(3.46) 

The  quantity  a is  the  constant 
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where  ax  ia  the  distance  between  mesh  points. 

Evidence  that  the  random  phase  screens/  generated  in 
the  manner  described  above/  have  the  desired  statistical 
properties  is  provided  by  the  phase  structure  function  results 
shown  in  Fig.  3.9.  The  theoretical  value  for  the  phase 
structure  function  for  the  spectrum  given  in  Eg.  (3.38)  is 


Vr)  “ " *<*2})2) 


2.61  k2  cj  Az  K“5/3  .6 
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<K0  r/2)V6  K5/6(K0  r) 


(3.48) 


which  in  the  limit  of  an  infinite  outer  scale  reduces  to 

( r ) - 2.91  k2  C2  A2  r5/3  ® (3.50) 

The  dashed  curved  line  and  the  straight  line  in  Fig.  3,9  are 
phase  structure  function  results  calculated  from  Eqs.  (3.48) 
and  (3.50).  The  somewhat  wicgly  solid  line  is  the  phase 
structure  function  result  obtained  by  averaging  over  100 
realizations  of  phase  screen  samples  generated  by  a fast 
Fourier  transform  evaluation  of  the  relation  given  in  Eq. 
(3.39).  In  these  runs  we  used  a 64  x 64  mesh  with  1.106  cm 
between  mesh  points  and  an  outer  scale  length  of  two  meters. 

3.4.6  Medium  Storage  and  Update  — At  each  (,  step 
in  the  propagation  calculation  we  store  the  absorptive  contri- 
bution to  P (although  this  is  done  in  large  core  memory  on  the 
CDC  7600  computer  that  we  are  using,  it  could  also  be  done  by 
writing  or.  the  disk).  The  effect  of  a transverse  wind  and 
of  slewing  is  taken  into  account  in  this  storage  process  by 
first  performing  a mesh  translation  on  the  new  and  old  data. 
The  data  are  moved  a distance  equal  to  the  distance  that  the 
wind  and  slewing  moves  the  heated  air  between  pulses.  Hence, 
the  stored  values  of  l'abB  are  those  that  will  be  seen  by  the 
next  pulse. 

Rather  than  calculate  l turb  anew  for  each  pulse  in  a 
multipulse  sequence,  we  calculate  a single  set  of  phase 
screens  at  the  outset  and  store  them  ft?  later  use.  The 
effect  of  the  wind  and  slewing  is  taken  into  account  by  trans- 
lating these  phase  screens  as  described  above. 

The  overall  logic  and  sequence  of  events  in  the  propa- 
gation simulation,  including  the  medium  Btorage  and  update 
processes  described  above,  are  indicated  in  the  flow  diagram 
given  in  Fig . 3.10. 
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SEPARATION  (r/l  105  cm  ) 


iy.  3.9.  Comparison  of  computer  calculated  and 
theoretical  phase  structure  function 
results. 
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AT  FOCAL  PLANE  7 


NO 


3 . 5 Target  Simulation 


The  essence  of  the  closed-loop  adaptive  optical  control 
concept  is  that  the  reflection  characteristics  of  typical 
targets  are  such  that  highlights  or  glints  will  exist  and 
that  the  return  from  these  regions  can  be  used  as  a reference 
to  sense  the  changes  in  the  transmitter  phase  distribution 
required  to  correct  for  atmospheric  phase  distortions.  Hence, 
we  approximate  the  target  by  a set  of  glints.  Each  glint  is 
represented  as  a gaussian  reflectivity  distribution,  i.e., 


M 

Reflectivity  Distribution  ■ S p <x)  (3.51) 
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m*l 


pm(x)  » a_  exp 
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(3.52) 


where  m is  the  number  of  glints  and  a . x , xm  are  the  strength, 
position  and  width  of  the  mcn  glint. 

The  initial  data  for  the  return  wave  are  determined  by 
multiplying  the  incident  laser  field  by  the  target  reflec- 
tivity function  given  in  Eg.  (3.51). 
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IV. 


SIMULATION  RESULTS 


4 . 1 Mirror  Simulation  Results 


The  major  portion  of  the  effort  on  the  mirror  simulation 
has  been  devoted  to  the  problem  of  determining  the  number  of 
mirror  actuators  required  to  yield  a good  fit  to  the  phase 
surfaces  that  must  be  generated  to  compensate  for  the  effects 
of  atmospheric  turbulence.  In  addition  to  the  computer  simu- 
lation results,  we  have  obtained  analytical  estimates  of  the 
effects  of  turbulence  on  peak  irradiance  at  a target,  both 
with  and  without  adaptive  compensation.  The  analytical  results 
obtained  for  the  effect  of  turbulence  on  peak  irradiance  without 
adaptive  compensation  are  discussed  in  Section  4.1.1.  Analy- 
tical and  computer  simulation  results  obtained  from  the  effect 
of  adaptive  compensation  with  a mirror  having  a finite  number 
of  degrees  of  freedom  (i.e.,  actuators)  are  discussed  in 
Section  4.1.2. 

We  also  have  done  some  work  on  the  dependence  of  mirror 
fit  on  the  number  of  mirror  actuators  for  the  case  of  phase 
conjugate  compensation  for  the  effects  of  thermal  blooming. 

These  results  are  discussed  in  Section  4.1.3. 

4.1.1  Effect  of  Turbulence  on  Peak  Irradiance  — An 
interpretation  of  theoretical  results  derived  for  the  ensemble- 
averaged  irradiance  distribution  of  a laser  beam  propagated 
through  atmospheric  turbulence5  indicatea  that  the  effect  of 
the  turbulence  can  be  accounted  for  by  introducing  an  appro- 
priately chosen  random  phase  screen  in  front  of  the  trans- 
mitter aperture.  Utilizing  thiB  result,  we  can  write  the 
Strehl  ratio  (SR)  associated  with  the  turbulence  in  the  form 
(recall  that  the  Strehl  ratio  due  to  a wave  distortion  is 
defined  as  the  ratio  of  the  peak  irradiance  in  the  presence  of 
the  distortion  to  the  peak  irradiance  of  the  undistorted  beam) 


Preceding  pege  blank 
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where  u^  is  the  aperture  field  distribution  and  D (x'-x")  is 

O (J)  — — 

the  phase  structure  function  of  the  random  phase  screen.  For 
Kolmogorov  turbulence  with  an  infinite  outer  scale,  the  appro- 
priate phase  structure  function  is 


(x  *-x") 
4,  - - 


2.91  k2 


(4.2) 


2 

where  is  the  index  structure  constant,  k is  the  wavenumber 
2 n / x and  z is  the  distance  to  the  target. 

For  a uniformly  illuminated  circular  aperture  of 
diameter  d,  the  Strehl  ratio  obtained  from  Eg.  (4.1)  is 
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with 
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where  ro  la  the  phase  cohoronce  length  defined  by  Fried 
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Numerical  reoults  for  the  Strehl,  ratio  defined  in  Eq.  (4.3) 
are  shown  in  Fig,  4.1  (curve  labeled  "no  adaptive  optics") 
for  values  of  [>,  up  to  15.  For  larger  values  of  [J, , the  follow 
ing  asymptotic  result  applies 


SR  % 4.407  r6/5  - 5.G92  if 9/5 , H • 15  (4.7) 
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Note  that  for  sufficiently  large  n , the  Strehl  ratio  is  jiven 

by  the  first  term  in  Fq.  (4.7)  and  from  Kq.  (4.5)  th- s c;.in  be 

2 

written  simply  as  the  ratio  (rQ/d)  , which  is  inversely  related 
to  the  number  of  phase  coherence  areas  contained  in  the  trans- 
mitter aperture. 

The  effect  of  having  an  adaptive  tilt  control  system 
that  dynamically  removes  the  tilt  component  of  the  turbulence 
distortion  can  be  accounted  for  by  replacing  the  structure 
function  D ^ in  Eq.  (4.1)  by  D(  - < (a  • (x'-x“))2>  where  a is 
a least  squares  estimate  of  the  instantaneous  tilt  vector. 
Utilizing  the  results  for  a derived  by  Fried,6  we  can  write 
the  following  expression  for  the  Strehl  ratio  obtained  for  a 
uniformly  illuminated  circular  aperture  having  adaptive  tilt 
control 


SR)with  tilt 
control 
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r.  (l-r.2) 


exp  -fi  f. 
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(4.8) 


Numerical  results  for  the  Strohl  ratio  defined  in  Eq.  (4.8) 

are  shown  in  Fig.  4.1  (curve  labeled  "with  tilt  control"). 

To  assist  in  the  use  of  Fig.  4.1  we  have  plotted  u as  a function 
2 

of  CN  z in  Fig.  4.2.  It  appears  that  tilt  control  alone  is  a 

reasonably  effective  compensation  technique  for  values  of  ('■ 

less  than  20.  At  10.6  urn  this  condition  will  be  satisfied 

in  many  scenarios,  wheroaa  at  3.8  1 1 m it  will  seldom  be 

satisfied  (e.g.,  we  note  from  Fig.  4.2  that  ('■  is  less  than  20 

2 - 10 

for  .7  m optics  at  10.6  urn  for  all  values  of  CN  z 2 x 10 

i/n  5 N — 15 

mi/  which,  for  a 4 km  path,  requires  that  z ••  2.5  x 10 

-2/3. 
cm  ) . 
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Fig.  4.2.  B versus  CN  for  homogeneous 
turbulence . 


4.1.2  Adaptive  Compensation  for  Turbulence  Effects  — 
Effect  of  Having  a Finite  Number  of  Actuators  — - It  was  noted 
in  the  previous  section  that  the  effect  of  atmospheric  turbu- 
lence on  the  ensemble-averaged  irradiance  distribution  can  be 
accounted  for  by  introducing  a random  phase  Bcreen  in  front 
of  the  transmitting  aperture.  The  appropriate  phase  distri- 
bution is  that  which  would  be  produced  by  a point  source 
located  on  the  target.  This  implies  that  an  ideal  adaptive 
optics  system  should  be  able  to  completely  compensate  for  the 

effects  of  the  turbulence  by  introducing  the  conjugate  phase 

+ 

distribution  as  a correction.  The  degree  to  which  this  is 
achieved  depends  on  whether  the  target  has  a point-glint 
scatterer  that  provides  the  requisite  phase  information  and 
on  the  degree  to  which  the  adaptive  optics  system  can  repro- 
duce the  desired  phase  distribution.  In  this  section  we 
shall  assume  that  the  target  has  a point-glint  scatterer  so 
that  the  ability  of  the  adaptive  optics  system  to  compensate 
for  the  turbulence  depends  solely  on  its  ability  to  reproduce 
the  desired  phase  distribution.  This  depends  largely  on  the 
number  of  degrees  of  freedom  possessed  by  the  correcting 
optics . 

To  a good  approximation,  we  can  represent  the  finite 
spatial  bandwidth  of  the  optics  by  writing  the  residual  phase 
error,  i.e.,  the  difference  between  the  desired  phase  distri- 
bution and  the  actual  distribution  obtained  from  the  optics, 
in  the  form 


ft 

In  this  discussion,  and  in  that  given  in  Section  4.1.1,  we 
assume  that  amplitude  scintillation  effects  arc  small.  If 
they  are  not  small,  the  effect  of  atmospheric  turbulence 
cannot  be  represented  simply  by  a phase  screen  but  rather 
we  must  also  introduce  a random  apodization. 

4* 

Subject  to  the  assumption  that  amplitude  effects  are  negli- 
gible . 


(4.9) 


(x)  = | (x) 


j dx  ' h (x  - x' ) <|j  (x ' ) 
j dx'  h(x') 


with 


h(x)  = 1,  j x j < b 

= 0,  | x | > b (4.10) 

In  this  representation  the  effect  of  the  finite  spatial  band- 
width of  the  optics  is  approximated  by  a spatial  averaging 
process,  i.e.  , we  assume  that  the  phase  correction  introduced 
by  the  adaptive  optics  at  a point  is  a uniformly  weighted 
average  of  the  desired  phase  correction  at  all  surrounding 
points  within  a circle  of  radius  b.  For  deformable  optics,  we 
can  relate  the  value  of  b to  the  number  of  mirror  actuators 

by  noting  that  the  spatial  spectrum  of  the  weighting  function 
h(x)  is  { 2 (Kb) /Kb] ^ , which  has  appreciable  magnitude  only  in 
the  domain  Kb  A-  3.  More  precisely,  we  have  found  that  good 
agreement  is  obtained  botween  the  theory  discussed  here  and  the 
mirror  simulation  results  discussed  later  if  Kb  « 2.72  where 
K is  the  cutoff  wavenumber  of  the  mirror.  (Note  that  the 
spectrum  of  the  weighting  function  h is  loss  than  .1  for  all 
Kb  ■ 2.72.)  Hence,  since.-  the  Nyquist  frequency  of  a mirror 

havjng  actuators  separated  by  a distance  v is  it /v  . we  let 

a a 


where  d is  the  mirror  diameter. 

The  Strehl  ratio  obtained  in  the  presence  of  the  residual 
phase  error  A*),  defined  in  Eq . (4.9)  is  given  by  Eq . (4.1)  with 

the  structure  function  D replaced  by 
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which 
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in  the  case  of  Kolmogorov  turbulence  with  an  infinite 
scale  can  be  written 


(4.13) 


with  u 


x - 


ilium!  n*.  A , r°  9 blr  Eq'  l4-2K  p«  « uniformly 

illuminated  oircular  aperture  of  diameter  d,  the  corre.pondino 

expression  for  the  strehl  ratio  obtained  with  adaptive  cor- 
rection  is 
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(4.15) 


We  have  evaluated  the  inteqr.l  in  Eq.  (4,1.)  numeri- 
cally and  have  found  that,  over  a wide  ran9e  of  value,  for  the 
parametc  a B and  b/d,  it  ia  approximated  extremely  well  by  the 
reeult  obtained  when  f(r,b/d)  i.  replaced  by  its  a.ymptotic 
result  for  b/d r,  1,  which  is 
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The  corresponding  expression  for  the  SR  is 

SR  'v,  exp  145  ft  ^ jj j ^ 

» exp^-,093  ft  n"5/6^  (4.17) 

wliere  we  have  used  the  relation  between  b/d  and  N given  in 

<9 

Eq.  (4.11).  The  discrepancy  between  the  approximation  in  Eq. 
(4.17)  and  the  more  exact  result  in  Eq.  (4.14)  is  significant 
only  for  very  large  values  of  ft  and  small  values  of  UQ.  This 
is  illustrated  by  Fig.  4.3  which  shows  the  results  obtained 
for  6 ■ 35.41  and  177,05.  For  values  of  N.  greater  than  5 the 

a 

discrepancy  is  less  than  one  percent. 

To  provide  a check  on  the  above  analytical  work,  we 
have  performed  two  types  of  mirror  simulations,  In  each,  the 
effect  of  a turbulent  phase  screen  placed  in  front  of  an 
aperture  was  compensated  by  a mirror  with  a finite  number  of 
actuators.  The  average  Strehl  ratio  obtained  from  a sequence 
of  five  independent  phase  screens  was  determined  as  a function 
of  the  number  of  actuators.  In  one  of  the  simulations  the 
deformable  mirror  software  described  in  Section  3. 3 was  used 
in  conjunction  with  a sinusoidal  multidither  COAT  algorithm  of 
the  type  used  previously  at  URL  in  our  COAT  servo  system 
studies.  In  the  other  simulation  the  mirror  was  modeled  by 
a segmented  mirror  with  piston  and  tilt  control  on  each  segment. 
The  piston  and  tilt  settings  for  each  segment  were  determined 
by  a least  squares  fit  to  the  random  phase  surface  over  the 
segment.  The  rosults  obtained  from  these  simulations  are 
summarized  in  Fig.  4.4.  The  agreement  between  the  piston  and 
tilt  resultB  and  the  theoretical  predictions  obtained  from 
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Fig.  4.3.  Comparison  of  Strehl  ratio  estimates. 


Eqs.  (4.14)  and  (4.17)  is  very  good  for  all  values  of  3 , 
which  gives  us  confidence  in  the  theoretical  results  given  in 

these  equations.  We  believe  that  these  results  are  repre- 
sentative of  those  that  would  bo  obtained  with  a deformable 
mirror  in  the  absence  of  the  2Nn  problem  discussed  below. 

The  results  actually  obtained  with  a deformable  mirror  will 
depend  on  the  degree  to  which  the  2Nti  problem  is  avoided, 

The  agreement  between  the  theory  and  the  deformable 
mirror-multidither  COAT  simulation  results  is  reasonably  good 
for  moderate  values  of  6 but  is  poor  for  large  values  of 
We  attribute  the  poor  results  obtained  at  large  values  of  3 to 
a 2Nir-type  of  behavior.  The  deformable  mirror  simulation  that 
was  used  in  these  runs  has  a 2Nit  correction  loop  that  introduces 
a 2 11  correction  whenever  the  phase  difference  between  actuators 
exceeds  four  radians.  The  intent  is  to  suppress  2 11  errors 
introduced  by  the  servo  system.  However,  if  the  phase  distortion 
that  is  to  be  removed  by  the  mirror  changes  by  more  than  four 
radians  between  actuators,  the  2 it  "correction"  introduced  by 
the  2Nit  correction  loop  is,  in  fact,  a 2 1 error  and  the 
mirror  performance  is  corresponding  degraded.  This  problem 
could  be  avoided  by  removing  the  2Nu  correction  loop  but  then 
we  would  be  faced  with  2Nn  servo  errors.  One  way  to  avoid 
this  problem  is  to  use  more  actuators  so  that  the  phase 
change  between  actuators  never  exceeds  four  radians.  Another 
way  would  bo  to  design  the  mirror  so  that  the  likelihood  of 
2 N is  errors  is  reduced,  in  which  case  the  2Nn  correction  loop 
could  be  eliminated.  More  work  clearly  needs  to  be  done  on 
this  problem. 

4.1.3  Compensation  for  Thermal  Blooming  Phase 
Prof  ilea  — As  a part  of  our  effort  to  provide  an  understand- 
ing of  the  reasons  for  the  relatively  meager  improvement  in 
beam  quality  realized  with  phase  conjugate  compensation,  we 
have  investigated  the  nature  of  the  phase  distribution  of  the 
return  wave  in  the  thermal  blooming  problem.  We  have  found 
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that  the  phase  of  the  return  wave  reflected  from  a qaussian 

glint  in  the  target  plane  does  not  vary  nearly  as  much  as 

one  would  expect  baaed  on  the  type  of  calculation  performed 

2 

by  Lincoln  Laboratory  in  their  optimization  studies. 

Specifically,  the  phase  correction  used  in  the  Lincoln 
Laboratory  work  is  of  the  form 

A*  - J sp  tn(l  + Nul)  / 1 dr'  |u(r.{,c2,0)  I2  (4.18) 

where  NQ  and  Nw  are  the  distortion  and  slewing  numbers  which, 
in  the  notation  used  in  this  report,  are  given  by  the  relations 
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yQ  is  the  polarizability  of  the  air,  r is  the  ratio  of  the 
specific  heats  at  constant  pressure  and  volume,  u is  the 
absorption  coefficient,  P is  the  average  power,  k is  the 
wavenumber,  f is  the  focal  distance,  v0  is  the  crosswind 
velocity,  P(j  is  the  ambient  atmospheric  pressure,  (>o  is  the 
e-folding  radius  of  the  laser  beam  (which  we  assume  to  have 
a qaussian  irradiance  distribution),  and  is  the  slew  rate. 
The  factor  t is  a constant  that  the  Lincoln  Laboratory  work 
indicates  is  nearly  equal  to  1.5  for  good  phase  compensation 
under  a variety  of  initial  beam  conditions. 

For  a gaussian  beam  truncated  outside  a circle  of 
diameter  3|,  , the  maximum  phase  excursion  predicted  by  Eq. 
(4.18)  is 
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Max.  Phase  Excursion  * A <}>  («•)  - A<|>  (-,u) 

Y ND  / X*  5 , 

- 5 fr  ln(l+NJ  / dr;  exp 

.856  Y Nn 

- « H ln(l+Nu)  (4.21) 

and  if  y = 1.5,  Nd  ■ 21.6,  and  Nw  * 4,  this  predicts  a maxi- 
mum phase  excursion  of  11.2  radians.  In  comparison,  we  have 
measured  a maximum  phase  excursion  for  the  field  reflected 
from  a gaussian  glint  under  equivalent  conditions  (for  ND  and 
NU))  to  be  only  about  two  radians.  We  believe  that  this  dis- 
crepancy is  caused  principally  by  the  extended  nature  of  the 
target  glint,  i.e.,  it  is  a gaussian  glint  rather  than  a point 
glintB,  and  that  this  discrepancy  supports  our  contention  that 
the  disappointingly  small  amount  of  correction  achieved  with 
phase  conjugate  compensation  is  related  to  the  isoplanatic 
problem  discussed  in  Section  2.1. 

4 . 2 Phase  Compensation  for  the  Effects  of  Turbulence 
and  'Thermal'  Blooming 

As  mentioned  earlier,  wo  have  obtained  phase  compen- 
sation data  only  for  the  case  iri  which  a phase  conjugate  control 
policy  is  used.  Since  all  return  wave  adaptive  optical  systems 
use  the  conjugate  phase  as  a correction  (see  the  discussion 
in  Section  2,1),  the  data  that  we  have  obtained  apply  to  all 
such  systems}  e.g.,  they  apply  to  conventional  phase  conjugate 
systems,  to  return-wave  multidither  systems,  and  to  all 
systems  based  on  a maximization  of  the  sharpness  of  an  image 
of  the  target. 

All  of  the  data  were  obtained  with  the  time-dependent 
propagation  code  described  in  Section  3.4.  Since  we  were 
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primarily  interested  in  the  correction  obtainable  under  cw 
conditions,  we  attempted  to  choose  pulse  repetition  frequen- 
cies that  yield  good  approximations  to  the  cw  blooming  condition. 
The  pulse  repetition  frequency  required  to  yield  conditions 
comparable  to  those  existing  with  cw  propagation  depends  on 
the  beam  diameter,  distance  to  the  focus,  wind  velocity,  and 
slew  rate.  Generally,  it  is  necessary  to  choose  a repetition 
rate  high  enough  that  an  air  parcel  moves  only  a fraction  of 
the  diffraction-limited  beam  diameter  between  pulses.  This 
condition  can  be  expressed  in  the  form 

( 1 * o;) 

Pulse  Repetition  Frequency  ■>  £ (4.22) 

where  v is  the  transverse  wind  velocity,  U is  the  slew  rate, 
f is  the  focal  distance,  k is  the  wavenumber  and  p is  a 
characteristic  beam  dimension  of  the  transmitted  beam  (e.g., 
the  e-folding  radius  of  a gaussian  beam). 

Figure  4.5  shows  a comparison  of  results  obtained 
with  the  time-dependent  code  and  a cw  code.  In  this  case 
we  plot  tho  Strohl  ratio  (peak  irradiance  in  the  presence  of 
distortion  divided  by  the  peak  irradiance  obtainable  without 
distortion)  of  an  untruncatod  gausBian  beam  focussed  at  2 km. 
Other  parameter  values  applicable  to  these  results  are  an 
averuge  power  level  of  50  kW,  a wavelength  of  10.6  urn,  an 
e-folding  radius  iQ  of  23  cm  (90%  of  the  energy  of  the  initial 
gauBsian  beam  contained  with  a 70  cm  diameter  circle),  an 
absorption  coefficient  of  2 x 10  ® cm"^,  a crosewind  velocity 
vo  of  10  m/sec,  and  a slew  rate  \i  equal  to  zero.  Apparently, 
a pulse  repetition  frequency  of  1000  Hz  or  more  is  required  to 
simulate  cw  propagation  under  those  conditions.  This  result 
is  in  qood  agreement  with  tho  condition  given  in  Eq . (4.22), 
which  predicts  that  vr  must  exceed  682  Hz.  As  an  aside,  it 
is  interesting  to  note  the  reduction  in  the  distorting  effect 
of  thermal  blooming  at  the  lower  repetition  frequencies.. 
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Fig.  4.5.  Strenl  ratio  versus  pulse  repetition  frequency  and  time 


It  becomes  increasingly  more  expensive  to  simulate  cw 
conditions  as  the  slew  rate  increases.  The  reason  for  this 
can  be  seen  from  the  data  in  Fig.  4.5.  Note  that  the  number 
of  pulses  required  to  reach  a steady  state,  be  it  the  cw 
state  or  not,  tends  to  increase  linearly  with  increasing  pulse 
repetition  frequency.  But  increasing  slew  rates  imply 
increasing  pulse  repetition  frequencies,  as  indicated  in  Eq. 
(4.22).  Hence,  we  must  propagate  more  pulses  to  reach  a 
steady  state  as  ft  increases.  In  our  calculations  we  have 
used  the  largest  pulse  repetition  frequency  consistent  with 
the  dual  requirements  that  we  simulate  the  cw  state  and  yet 
do  not  require  an  excessive  number  of  pulses.  Typically,  we 
have  used  no  more  than  20  pulses  for  a given  simulation.  At 
high  slew  rates,  this  has  resulted  in  some  departure  of  the 
multipulee  of  results  from  those  that  would  be  obtained  from 
a cw  code.  Principally,  the  difference  is  that  for  a given 
average  power  the  multipulse  results  tend  to  yield  Strehl 
ratios  that  are  30  to  401  higher  than  those  obtained  from  a 
cw  code.  The  same  type  of  power  optimization  curves  are 
obtained  except  that  the  peak  irradiance  and  the  transmitter 
power  required  to  yield  the  peak  are  higher  (see  Fig.  4.6 
for  example) . This  discrepancy  between  the  multipulse  and 
cw  results  is  not  important  in  our  studies  because  the  main 
goal  of  these  studies  is  to  ascertain  the  efficacy  of 
adaptive  optics  as  a technique  for  increasing  the  magnitude 
of  the  peak  irradiance  at  the  target,  not  to  establish  the 
precise  power  level  at  which  that  peak  occurs. 

In  order  to  determine  the  relative  amount  of  correction 
to  be  expected  for  turbulence  and  thermal  blooming  effects, 
we  have  made  runs  in  wnich  one  or  the  other  ot  these  effects 
was  absent.  in  addition,  we  have  obtained  data  for  adaptive 
compensation  of  the  combined  effects  of  turbulence  and 
thermal  blooming.  In  all  of  these  calculations  we  have  used 
a phase  conjugate  correction  algorithm  in  which  the  transmitter 
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phase  was  equated  on  a mesh-point-by-mesh-point  basis  to  the 
negative  of  the  phase  returned  from  the  target.  The  target 
was  represented  by  a single  gaussian  glint  having  an  e-folding 
width  equal  to  4 cm.  The  results  obtained  for  propagation  at 
a COj  laser  wavelength  of  10.6  pm  are  discussed  in  Section 
4.2.1  and  those  obtained  for  propagation  at  a DF  laser  wave- 
length of  3.8  im  are  discussed  in  Section  4.2.2.  In  these 
calculations  the  absorption  coefficient  was  assumed  equal  to 
2 10  0 cm  at  10.6  urn  and  4 x 10  cm  at  3.8  urn. 

4.2.1  Phase  Conjugate  Compensation  at  10.6  pm  — One 
of  the  more  interesting  results  from  these  studies  iB  the 
observation  that  much  better  thermal  blooming  compensation 
is  obtained  with  untruncated  gaussian  beams  than  with  trun- 
cated gaussian  beams.  An  example  of  the  results  obtained  with 
an  untruncated  gaussian  beam  is  shown  in  Fig.  4.7,  which 
applies  to  the  propagation  of  a beam  having  an  e-folding  radius 
of  23  cm  (90%  of  the  energy  in  the  initial  irradiance  distri- 
bution is  contained  within  a circle  having  a 70  cm  diameter) , 

a slew  rate  of  20  mrad/sec,  a focal  distance  of  2 km,  and 

2 

an  index  structure  constant  CN  » 0 (i.e.,  no  turbulence). 

Phase  conjugate  compensation  increases  the  peak  irradiance 
by  a factor  of  2.3  and  the  transmitter  power  at  which  the 
peak  occurs  is  increased  by  a factor  of  two.  For  comparison, 
the  results  obtained  for  a gaussian  beam  truncated  at  a beam 
radius  of  35  cm,  with  all  other  beam  and  scenario  parameters 
identical  to  those  discussed  above,  are  shown  in  Fig.  4.8. 

In  this  case,  phase  conjugate  compensation  increases  the 
peak  irradiance  only  by  about  30%.  We  believe  that  the 
difference  in  phase  compensation  achieved  with  and  without 
beam  truncation  can  be  explained  by  the  isoplanatic  problem 
discussed  in,  Section  2.1.  As  Bhown  in  Fig.  4.9,  the  fields 
within  the  circle  of  beam  truncation  returned  from  points 
within  the  domain  of  an  extended  glint,  such  as  the  gaussian 
glint  used  in  our  numerical  studies,  pass  through  air  that 
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Pig.  4.B.  Phase  conjugate  compensation  of  a 
truncated  gaussian  beam. 


Fig.  4.9.  Isoplanatic  problem. 


has  much  larger  density  gradients  than  the  fields  returned 
outside  of  the  circle  of  beam  truncation.  Hence,  the 
distorting  effect  that  a finite  glint  has  on  the  phase  of  the 

return  wave  Is  less  severe  outside  of  the  region  of  truncation 
than  it  is  within  this  region  and,  therefore,  better  compen- 
sation should  be  obtained  with  infinite  beams  than  with 
truncated.  In  fact,  it  appears  that  for  the  infinite  beam 
case  shown  in  Fig.  4.7  almost  all  of  the  factor  of  2.3 
increase  in  target  irradiance  is  obtained  from  fields  outside 
of  a 70  cm  diameter  circle. 

Results  for  phase  conjugate  compensation  of  the 
combined  effects  of  turbulence  and  thermal  blooming  are  shown 
in  Figs.  4.10,  4.11  and  4.12.  All  of  these  results  apply 
to  truncated  gaussian  beams  propagated  through  a single 
realization  of  a random  medium;  i.e.,  no  ensemble-averaging 
over  different  random  medium  realizations  was  done.  Generally, 
these  results,  and  others  that  we  have  obtained,  indicate 
that  phase  conjugate  compensation  yields  good  compensation 
for  the  effects  of  turbulence  but  very  little  for  the  effects 
of  thermal  blooming.  Fortunately,  the  ability  t.o  compensate 
for  the  effects  of  turbulence  is  apparently  not  degraded  by 
the  presence  of  thermal  blooming  effects  even  though  phase 
conjugate  compensation  iB  not  very  effective  as  a thermal 
blooming  compensation  technique.  The  not  improvement  in  beam 
quality  obtained  with  phase  conjugate  compensation  is  most 
impressive  for  the  conditions  depicted  in  Fig.  4.12  which 

applies  to  a 1,2  in  beam  propagated  through  moderately  strong 

- 1 5 

turbulence  having  an  index  structure  constant  of  6 x 10 
-2/3 

cm  . This  is  to  bo  expected,  of  course,  since  turbulence 
effects  are  most  pronounced  in  the  case  of  largo  beams  and 
high  turbulence  levels. 

The  severity  of  the  thermal  blooming  problem  for  cw 
propagation  at  10.6  urn  has  engendered  a considerable  amount, 
of  interest  recently  in  multi  pulse  propagation  in  which  the 
time  between  pulses  is  sufficient,  to  allow  the  local  cross- 
wind  velocity  to  convect  partially  the  heated  air  from  the 
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region  occupied  by  the  beam.  This  is  a favorable  situation 

tor  thermal  blooming  compensation  because  the  blooming  under 
these  conditions  tends  to  occur  closer  to  the  transmitter 
and  is,  therefore,  more  nearly  a pure  phase  effect  which  can 
be  corrected  by  a compensating  phase  distribution  on  the 
transmitted  wave.  Moreover,  the  deleterious  effects  of  the 
isoplanatic  problem  discussed  in  Section  2.1  are  greatly 
alleviated  so  that  we  should  obtain  significantly  better 
phase  compensation  results  in  the  multipulse  mode  than  we  have 
for  cw  propagation.  This  is  illustrated  by  the  results  shown 
in  Pig.  4.13.  TheBe  results  apply  to  phase  conjugate  compen- 
sation for  the  residual  effects  of  thermal  blooming  on  the 
propagation  of  12  kJ  pulses  of  10.6  um  energy  radiated  in  a 
truncated  gaussian  beam  having  a diameter  of  70  cm  and  an 
e-foldinq  radius  of  23  cm  (90%  of  the  energy  in  the  gaussian 
contained  within  a circle  having  a 70  cm  diameter).  The 
crosswind  velocity  is  10  m/sec  and  the  beam  is  being  slewed 
at  20  mrad/soc.  The  results  are  plotted  as  a function  of 
average  power,  which  is  related  to  the  pulse  repetition  rate 
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(4.23) 


The  effect  of  thermal  blooming  within  a pulse  is  nog  looted 

in  this  calculation.  This  should  be  relatively  small  for 

short  pulse  lengths.  For  the  conditions  pertaining  to  the 

results  in  Fig.  4.13,  the  minimum  pulse  length  consistent  with 

having  a peak  irradiance  below  the  sea  level  air  breakdown 
7 2 

threshold  of  10  watts/cm  is  10  hboc.  (TIub  assumes  that 
the  Strehl  ratio  is  not  degraded  by  thermal  blooming  effects.) 
The  significant  increase  in  target  irradiance  obtained  with 
phase  compensation  verifies  that  multi  pulse  propagation  is 
a very  favorable  mode  of  operation  for  phase  compensation . 
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4.2.2  Phaso  Conj  ligate  Compensat  ion  a t_ 3 . 8_  j_i_m  — The 
phase  conjugate  correction  results  obtained  at  3.6  urn  are 
consistent  with  those  obtained  at  10.6  pm  in  that  wo  have 
observed  that  significant  correction  is  achieved  for  turbu- 
lence effects  but  that  very  little  correction  is  realized 
for  thermal  blooming.  The  impressive  amount  of  correction 
achievable  for  turbulence  effects  in  the  absence  of  thermal 
blooming  is  illustrated  by  the  Strahl  ratio  results  given 
in  Figs.  4.14  and  4.15  which  apply  to  a case  in  which  the 
index  structure  constant  is  cm”2^3.  The  results  in 

Fig.  4.14  pertain  to  a gaussian  beam  with  an  e-folding  radius 
of  13.8  cm  (90%  of  the  energy  in  the  initial  irradiance 
distribution  is  contained  within  a circle  having  a 40  cm 
diameter)  truncated  outside  a 70  cm  diameter  circle.  The 
path  length  is  2 km  and  the  slow  rate  is  20  mrad/sec.  The 
results  in  Fig.  4.15  apply  to  a gaussian  beam  with  an  e-folding 
radius  of  40  cm  (90%  of  the  energy  in  the  initial  irradiance 
distribution  is  contained  within  a circle  having  a 1.2  m 
diameter)  truncated  outBide  a 1.2  m circle.  The  path  length 
in  thiB  case  is  4 km  and  the  slew  rate  is  5 mrad/sec.  In 
both  cases  the  Strahl  radio  obtained  with  correction  eventually 
attains  a value  within  10  to  20%  of  the  freo-space  value. 
Although  the  results  in  Fig.  4.15  exhibit  a rather  slow 
increase  in  Strehl  ratio  with  successive  pulses  until  the 
twelfth  pulse,  we  believe  that  an  actual  phase  conjugate 
system  would  increase  the  Strehl  ratio  much  faster.  The 
Blow  build-up  seen  in  Fig.  4.15  is  probably  caused  by  the 
fact  that  we  did  not  iterate  the  correction  process  at  each 
step*  rather,  a single  phase  measurement  was  made  based  on 
the  return  of  a single  reference  pulse  propagated  just  prior 
to  the  high  power  pulse.  The  single  pulse  leferenco  is  adequate 
when  the  beam  is  in  a nearly  converged  state,  as  it  is  later 
on  in  the  pulse  train,  but  does  not  give  an  accurate  phase 
estimate  early  in  the  pulse  sequence.  Early  in  the  pulse 
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train  the  return  comes  from  a more  extended  portion  of  the 

target  because  of  the  broadness  of  the  transmitted  beam  and 
the  very  significant  fluctuation  in  the  amplitude  of  that 
beam  in  the  target  plane. 

Results  for  phase  conjugate  correction  of  the  combined 
effects  of  turbulence  and  thermal  blooming  are  shown  in  Figs. 

4.16  and  4.17.  The  results  in  Fig.  4.16  apply  to  a truncated 
gaussian  beam  with  an  e-foiding  radius  of  23  cm  and  a trun- 
cation circle  diameter  of  70  cm.  This  beam  is  being  slewed  at 
20  mrad/sec.  The  results  in  Fig.  4.17  portain  to  a truncated 
gaussian  beam  with  an  e-folding  radius  of  40  cm  and  a trun- 
cation circle  diameter  of  1.2  m (same  degree  of  truncation  as 
the  70  cm  beam  used  in  Fiq.  21).  This  beam  is  beinq  slewed 
at  only  5 mrad/sec.  In  addition  to  demonstrating  the  ability 
of  phase  conjugate  adaptive  systems  to  achieve  substantial 
corrections  for  the  effects  of  turbulence#  these  results 
illustrate  the  advantage  of  using  large  beams  in  thermal 
blooming  situations.  Even  though  the  slow  rate  in  Fig.  4,17 
is  four  times  slower  than  that  in  Fig.  4,16#  the  peak 
irradianco  obtained  without  correction  is  increased  by  a 
factor  of  five.  Likewise#  the  peak  irradiance  obtained  in 
Fig.  4.17  with  phase  conjuqate  correction  is  incroased  by  a 
factor  of  three  beyond  the  comparable  result  obtained  in 

Fig.  4.16.  This  improvement  in  peak  irradiance  can  bo 
attributed  to  the  ineroaard  size  of  the  beam  used  in  Fig. 

4.17  (1,2  in  optics  versus  7C  cm  optics). 
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Fig.  4.16.  PhaBe  conjugate  compensation  at  3.8  urn, 
70  cm  beam. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

5 . 1 Conclusions 

The  most  important  conclusion  that  we  have  reached  on 
the  basis  of  the  work  done  on  this  contract  is  that  adaptive 
optical  systems  employing  any  of  the  various  types  of  return- 
wave  control  algorithms  are  not  likely  to  yield  very  singi- 
f leant  compensation  for  the  effects  of  cw  thermal  blooming 
unless  the  target  scattering  occurs  over  regions  having 
extent  consideiably  smaller  than  the  diffraction-limited  beam 
size.  As  a matter  of  fact,  some  of  the  data  we  have  obtained 
indicate  that  under  some  conditions  such  adaptive  systems 
may  actually  decrease  target  irradiance  instead  of  increasing 
it.  The  analysis  in  Section  2.1  of  this  report  indicates  that 
the  problem  with  the  return-wave  control  algorithm  originates 
from  the  fact  that  it  inherently  tends  to  use  the  conjugate 
of  the  phase  of  the  field  returned  from  the  target  as  the 
phase  correction.  This  phase  conjugate  control  policy  works 
well  as  long  aB  the  target  acts  as  a point  scatterer  or  a 
discrete  collection  of  such  scatterers.  However,  if  the 
scattering  originates  from  extended  scatterers  rather  than 
point  scatterers,  there  is  an  isoplanatic- type  of  effect 
wherein  the  return  phase  becomes  garbled  by  the  interference 
between  returns  from  the  various  parts  of  the  scatterer. 

This  iBoplanatic  effect  is  of  no  consequence  for  moat  turbu- 
lence compensation  situations  since  the  isoplanatic  area  is 
usually  larger  than  the  principal  scattering  region.  Hence, 
return-wave  systems  yield  good  compensation  for  turbulence 
effects.  For  cw  thermal  blooming,  however,  our  computer 
results  indicate  that  the  isoplanatic  region  is  smaller  than 
the  diffraction-limited  spot  size.  We  baso  this  conclusion 
on  the  fact  that  very  little  compensation  was  obtained  for 
cw  thermal  blooming  effects  with  glints  having  an  extent 
approximately  equal  to  the  diffraction-limited  spot  size. 
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It  should  be  noted  that  the  above  conclusions  apply 
to  any  adaptive  optical  system  employing  a return  wave-type 
of  control  algorithm,  not  ju3t  the  explicit  phase  conjugate- 
type  of  system.  As  discussed  in  Section  11,  this  includes 
systems  that  operate  by  maximizing  some  functional  of  the 
irradiance  distribution  obtained  by  imaging  the  target;  e.g., 
maximization  of  the  energy  through  a pinhole  for  a point 
glint  target  or  maximizing  the  "sharpness"  of  an  image  of  an 
extended  target. 

With  regard  to  turbulence  compensation,  we  have  found 
that  return-wave  adaptive  optical  systems  yield  significant 
corrections.  Typically,  it  is  possible  to  obtain  almost 
diffraction-limited  target  intensities  even  for  relatively 
strong  turbulence.  The  amount  of  improvement  is  particularly 
impressive  at  3.8  p m where  the  turbulence  degradation  is  quite 
severe  without  adaptive  compensation;  e.g.,  we  have  obtained 
a factor  of  seven  increase  in  peak  target  irradiance  for  a 
70  cm  diameter,  3.8  |im  beam  propagated  over  a 4 km  path  having 
a turbulence  level  equivalent  to  an  index  structure  constant 
of  cf.  * 10  cm  ' (moderate  turbulence  on  a low  altitude 

N 

path) . 

In  the  results  discussed  above,  and  in  all  the  other 
results  in  which  the  full  propagation  simulation  was  used, 
we  hiive  not  simulated  the  effect  of  havinq  fewer  mirror 
actuators  than  mesh  points.  We  have  simply  equated  the  phase 
correction  to  the  conjugate  of  the  return  wave  on  a niesh- 
point-by-mesh-point  basis.  However,  we  have  obtained  quanti- 
tative estimates  of  the  degrading  effect  of  having  fewer 
actuators  than  meBh  points  from  both  theoretical  analysis 
and  a simplified  computer  simulation  in  which  the  effect  of 
turbulence  was  represented  by  a single  phase  screen  placed 
in  front  of  the  transmitter  (see  Section  4.2  for  details). 
Using  the  theory  developed  in  Section  4.2.2,  which  yields 
results  that  agree  well  with  the  computer  simulation  results, 
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we  obtain  the  following  condition  on  the  number  of  actuators 
per  unit  area  of  mirror  surface  required  to  yield  a Strehl 
ratio  of  .9  or  greater  (i.e,,  the  degradation  in  Strehl  ratio 
due  to  the  finite  number  of  mirror  actuators  is  < .10) 

Number  of  Actuators  , , , 

Required  per  Unit  L .529(k*  cj  z)0/3  (5.1) 

Area 

We  estimate,  for  example,  that  approximately  .037 
actuators  per  cm  of  primary  mirror  area  are  required  to 
yield  results  comparable  to  the  mesh-point-by-mesh-point 
results  for  the  3.8  urn  case  cited  above,  i.e,,  approximately 
143  actuators  are  required  for  70  cm  optics. 

In  the  case  of  multipulse  propagation,  we  have  found 
that  adaptive  optical  systems  employing  return  wave-type  of 
control  algorithms  can  significantly  increase  the  target 
irradiance  in  thermal  blooming  situations.  The  amount  of 
correction  attained  varies  depending  on  the  conditions  but 
can  approach  a factor  of  three  in  some  cases.  The  reason  for 
the  significantly  better  correction  obtained  for  multiple 
pulse  thermal  blooming  effects  is  principally  that  the  distort- 
ing effect  of  the  thermal  blooming  tends  to  move  back  away 
from  the  target  and  closer  1.0  the  transmitter.  This  has  the 
effect  of  increasing  the  isoplanatio  area  at  the  target,  and 
thus  alleviates  the  problem  that  occurs  with  cw  thermal 
blooming . 

Although  we  have  not  yet  obtained  computer  simulation 
results  for  an  adaptive'  optical  system  employing  an  outgoing- 
wave  multidither  control  algorithm,  we  have  developed  an 
approach  to  the  problem.  This  approach  is  a modified  version 
of  what  has  been  done  at  Lincoln  Laboratory  by  Herrmann  and 
Bradley.  Since  an  outgoing-wave  multidither  control  algorithm 
is  basically  a max  i m i /.alt  on  process,  it  can  bo  represented 
by  a function  maximization  routine.  In  contrast  to  the 
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Lincoln  Laboratory  work  in  which  the  target  irradiance  was 
directly  maximized  with  respect  to  the  phase  correction,  we 

recommend  a simulation  which  more  closely  models  the  operation 
of  an  actual  adaptive  optical  system.  To  accomplish  this,  we 
would  use  a time-dependent  propagation  code  and  would  maximize 
the  integrated  irradiance  of  the  field  returned  to  the  trans- 
mitter from  the  target.  There  are  two  principal  differences 
between  this  approach  and  that  taken  by  Lincoln  Laboratory. 
First,  we  would  maximize  the  integrated  irradiance  of  the  field 
returned  from  the  target  rather  than  the  target  irradiance 
iirectly.  This  will  allow  us  to  determine  the  sensitivity  of 
the  maximization  proceBB  to  the  size  of  the  receiver  and  the 
nature  of  the  target  glint.  Secondly,  use  of  the  time- 
dependent  code  will  allow  us  to  determine  whether  the 
instantaneous  maximum  provided  by  an  outgoing-wave  multi- 
dither  algorithm  yieldB  a long-term  global  maximum  irradiance 
at  the  target  or  instead  simply  yields  a secondary  maximum. 

Put  in  another  way,  we  would  be  able  to  determine  whether 
it  is  possible  to  end  up  at  a local  (secondary)  maximum 
rather  than  the  global  maximum  when  we  instantaneously  maxi- 
mize the  integrated  irradiance  of  the  return  fieldj  i.e., 
could  we  have  reached  a higher  maximum  had  we  taken  a 
different  path  in  time  (via  a different  control  policy)? 

5 . 2 Recommendations 

Dividing  the  subjects  again  on  the  basts  <«t  the  type 
of  control  algorithm  used,  wo  recommend  the  following  addi- 
tional work  on  the  adaptive  optics  simulation  problem: 

Return  Wave  Systems  — Generally,  moi e data  are  needed 
for  a broader  set  of  system,  atmospheric,  and  target  para- 
meters than  have  been  investigated  to  this  date.  These 
parameters  include  the  mirror  diameter,  type  of  mirror 
(segmented  or  deformable),  number  of  actuators,  initial 
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Lriadiancu  distribution  (including  distributions  typical  of 
actual,  high  power  lasers),  target  range  and  slew  rate,  and 

target  glint  size.  Of  particular  importance  in  the  cw 
thermal  blooming  problem  is  the  investigation  of  the  degree 
of  compensation  achievable  as  a function  of  target  glint  size. 
Those  results  should  be  correlated  with  known  properties  of 
real  targets  in  order  to  arrive  at  an  authoritative  assess- 
ment of  the  offects  of  target  glint  size  on  return-wave 
systems , 

Outgoing-Wave  Multidither  Systems  ■—  As  indicated 
earlier,  we  recommend  that  Buch  systems  be  simulated  by  a 
function  maximization  routine  in  conjunction  with  a time- 
dependent  propagation  code.  Initially,  at  least,  the  main 
intent  of  such  a simulation  should  be  to  establish  whether 
an  outgoing-wave  multidither  system  yields  a truly  global 
maximization  for  the  irradiance  at  the  target  or  instead 
yields  simply  a secondary  maximum.  The  effect  of  target 
glint  size  and  receiver  size  on  the  maximization  process  will 
be  of  particular  interest  in  these  studies.  Other  parameters 
that  should  be  varied  are  mirror  diameter,  type  of  mirror, 
number  of  actuators,  initial  irradiance  distribution,  and 
target  range  and  slow  rate. 
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APPENDIX  A 


RECIPROCITY  CONSIDERATIONS 


The  equivalence  of  the  transmitting  and  receiving 
patterns  of  an  antenna  is  a consequence  of  a general  reci- 
procity relation  satisfied  by  the  solution  of  the  Maxwell 
equations.  For  the  propagation  problems  of  concern  in  this 
report,  the  solution  of  the  Maxwell  equations  can  be  derived 
from  the  solution  of  a reduced  wave  equation.  The  desired 
reciprocity  relation  is  obtained  by  considering  the  case  where 
there  is  a delta  source  located  at  one  of  two  points  x^ 

(i  1,2)  within  a volume  V bounded  by  u surface  S 

|v2  + k2  n 2 ( x ) ] ui(x)  - Mx  - xi)  (A.l) 

We  assume  that  the  refractive  index  distribution  within  V 
and  the  boundary  conditions  on  S are  independent  of  the  location 
of  the  source.  If  we  multiply  the  equation  for  (x)  by  Uj  (x) 
and  that  for  u2<x)  by  Uj (x) , subtract  the  results,  and  inte- 
grate over  V,  wo  obtain 


U1  <V 


(A.  2) 


The  surface  integral  is  obtained  by  applying  Green's  theorem 

2 2 

to  the  volume  integral  ot  V Uj  - Uj  v u2< 

For  a wide  range  of  problems,  including  those  of 
interest  in  this  report,  we  can  approximate  the  boundary 
conditions  on  S by  the  impedance  boundary  condi tion 


b 


m 


- 0, 


x on  S 


(A.  3) 


where  u is  a constant  that  depends  on  the  surface  impedance 

Preceding  page  blank 


of  the  boundary.  Hence,  the  surface  integral  in  Eq.  (A. 2) 
is  zero  and  we  obtain  the  following  reciprocity  relation 

(Xj ) **  (A.  4) 

This  relation  states  that  the  field  produced  at  x2  by  n source 
at  x^  is  the  same  as  the  field  produced  at  x^  by  a source  at 
x^.  We  note  that  the  derivation  of  Eq.  (A. 4)  did  not  require 
any  assumptions  regarding  the  nature  of  the  refractive  index 
function  n(x)  in  Eq.  (A.l).  Hence,  the  reciprocity  relation 
we  have  derived  is  valid  regardless  of  the  state  of  homo- 
geneity of  the  medium  between  s^  and  x2> 

To  prove  the  equivalence  of  the  transmitting  and 
receiving  patterns  of  an  antenna  immersed  in  an  inhomoqeneous 
medium,  we  need  now  merely  consider  the  situation  shown  in 
Fig.  A.l.  The  antenna  is  depicted  here  as  a lens.  According 
to  the  reciprocity  relation  given  in  Eq.  (A. 4),  the  field 
radiated  to  the  point  x2  from  a point  source  at  x^  behind  the 
lens  is  equal  to  the  field  received  at  x^  from  a point  source 
at  x2*  This  implies  that  the  transmi  '.ng  and  receiving 
patterns  of  the  lens  are  the  same  because,  as  shown  in  Fig. 
A.l,  the  field  at  in  the  receiving  case  bears  the  same 
geometrical  relationship  to  the  center  of  the  receiving 
pattern  (indicated  by  the  dashed  lines  in  Fig,  A.l),  as  does 
the  field  at  x2  to  the  center  of  the  transmitted  pattern 
(indicated  by  the  double-dashed  line  in  Fig.  A.l). 

As  a corollary  to  the  equivalence  between  the  trans- 
mitting and  receiving  patterns  of  an  antenna,  we  note  that  if 
one  images  a point  source  located  at  x2,  the  point  of  maximum 
irradiance  in  the  image  plane  is  the  ideal  place  to  locate 
a point  source  if  it  is  desired  to  deliver  maximum  power  to 
x2  on  transmit.  To  prove  this,  we  recall  from  Eq.  (A. 4) 
that  the  field  at  x2  radiated  by  a source  at  x^  is  the  samo 
as  that  which  occurs  when  the  source  and  observation  points 
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arc  interchanged.  Hence,  the  fact  that  we  locate  the  point 
Bource  at  the  point  of  maximum  irradiance  from  the  source  at 
X£  implies  that  the  field  radiated  to  £j  will  be  maximized 
because  the  transmitting  and  receiving  patterns  of  the 
antenna  are  equivalent  and  a point  source  located  at  the 
maximum  of  one  pattern  will  likewise  produce  at  maximum  at 
the  reciprocal  point.  Basically*  we  are  saying  that  maximum 
power  is  delivered  to  x2  when  the  transmitter  is  pointed  at  the 
apparent  position  of  Xj . This  result  has  wide  application  in 
communication  and  radar  systems  where  it  is  often  necessary 
to  deliver  energy  to  a point  specified  by  an  imaging  system. 

In  addition  to  the  equivalence  between  the  transmitting 
and  receiving  patterns  of  an  antenna  and  the  pointing  corol- 
lary discussed  above,  the  reciprocity  relation  given  in  Eq. 

(A. 4)  has  the  obvious  consequence  that  the  fluctuations  of  a 
spherical  wave  propagating  in  a random  medium  are  reciprocal 
regardless  of  whether  the  random  inhomogenej ties  are  spatially 
homogeneous.  By  reciprocal  wc  mean  that  the  fluctuations  at 
a point  £2  in  the  field  radiated  from  a spherical  wave  source 
at  a point  aro  the  same  as  those  observed  at  x^  when  the 
spherical,  source  is  at  £2*  This  means,  for  instance,  that 
oarthbound  measurements  of  the  fluctuations  of  a spherical  wave 
emanating  1'rom  a point  Boutco  outside  the  earth's  atmosphere 
can  be  used  to  predict  the  fluctuations  that,  would  occur  In 
the  reciprocal  situation*  namely,  an  ear t.hbound  point  aourco 
radiating  to  the  same  point  in  space. 

Unfortunately,  the  reciprocal  relation  between  the 
fluctuations  of  the  spherical  wave  fiolds  in  random  media 
generally  does  not  apply  to  the  fluctuations  of  fields  radiated 
by  more  complicated  source.  For  oxample,  consider  the  source- 
medium  configuration  shown  in  Figs.  A.*a  and  b.  The  source 
in  this  figure  is  composed  of  a point  aourco,  lens  combination, 
which  is  gonorically  equivalent  to  most  laser  source.  With 
this  type  of  Bource  we  change  the  location  of  the  ions  when 
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the  source  and  observation  points  are  interchanged.  In 

effect,  this  means  that  the  refractive  index  function  n(x) 
in  Eq.  (A.l)  is  not  the  same  when  the  ''source"  is  at  x^, 
as  it  is  when  the  "source"  is  at  Xj.  Consequently,  the  terms 
involving  the  refractive  index  do  not  completely  cancel  when 
we  follow  the  procedure  outlined  earlier  for  the  derivation 
of  the  reciprocity  relation,  given  in  Eq.  (A. 4).  There  is  an 
additional  term  in  Eqs.  (A. 2)  and  (A. 4)  that  destroys  the 
simple  reciprocal  relationship  that  exists  for  point  source 
fields. 
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